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A stereotactic device 
for use with an image, which 
permits an imaging scan to 
be taken with reference to a 
personal coordinate system 
(PCS) that is independent of 
a machine coordinate system 
(MCS). Methods using the 
device to obtain imaging 
scans (28) arc described such 
that the imaging scans are 
superimposing even if taken 
at different time periods using 
the same or a different image. 
The device comprises a frame 
(10) that can be rcproducibly 
positioned on a subject, 
and which is equipped with 
non- invasive affixing means 
(12, 14) and localizing means 
(16A, 16B, 18A, 18B. 20A, 
20B) that provide the PCS. 
The device and method steps 
are particularly well suited 
for routine initial or follow-up 
examinations, pre-surgical 
planning and post-surgical 
evaluation. 
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VERSATILE STEREOTACTIC DEVICE AND METHODS OF USE 

5 

Cross Reference to Related Application 
This application claims the priority date of Provisional 
U.S. Application Serial No. 60/002,497, filed August 18, 1995, 
the disclosure of which is incorporated in its entirety by 
10 reference herein. 

1 . Field of the Invention 

The invention pertains to a versatile stereotactic device 
useful in a number of methods, including numerous modes of 

15 medical imaging. More particularly, the device and methods of 
the present invention relate to a non- invasive stereotactic 
method of reproducibly imaging portions of a patient's body, 
such as the patient's head and portions of the patient's spine 
in the proximity of the head. Thus, imaging modalities, 

20 including magnetic resonance (MR) imaging, magnetic resonance 
spectroscopy, computer-aided tomography (CT) , positron emission 
tomography (PET) , single photon emission computed tomography 
(SPECT) , electroencephalography (EEG) or magnet ©encephalography 
(MEG) and the like can be used to monitor, diagnose, or detect 

25 pathologic conditions and to follow their development, 
progress, arrest, or remission. The device and methods of the 
invention are especially applicable to permitting more routine 
head examinations, pre-surgical planning and providing post- 
surgical evaluations and prognoses. 

30 

2 . Background of the Invention 

In head examinations involving magnetic resonance imaging, 
computer-aided tomography and other such techniques, it is 
desirable to have a well-defined, reproducible coordinate 
35 system to record and/or compare the locations and sizes of 
lesions, tumors and other structures of interest. Though there 
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are a number of known devices and techniques for potential 
application to these types of examinations, these known devices 
and techniques are generally not suitable for ^routine" office 
examinations, in which factors such as ease of use, speed of 
5 use, comfort, cost, accuracy and reproducibility are of major 
consideration. Indeed, existing devices are often heavy, 
unwieldy, cumbersome and require that the devices be affixed to 
the subject using pins, screws, bolts, brackets, staples and 
the like. 

10 Several methods have been proposed to find the relative 

position of a scan "slice" by. using anatomical landmarks. In 
these methods, the size and position of predetermined 
anatomical structures, such as the lateral end of the internal 
auditory canal, are used as reference points to help locate and 

15 compare lesions and other features of interest. See, e.g., 
Tan, K.K. et al., in J. Neurosurg. (1993) 79:296-303. A 
problem with this technique, however, is that the image 
resolution in the scan slice direction (i.e., the z-direction) 
is poor compared with the resolution in the scan in-plane (x-y 

20 direction) . Because of this poor resolution, it is difficult 
to make precise positional determinations. Moreover, these 
methods also require a degree of anatomical knowledge which may 
strain the capabilities of the average MR technologist. 

As an alternative, the art has developed devices, such 

25 as frames and w halos," to facilitate positioning for 
stereotactic surgery. These devices are rigidly affixed to the 
patient being imaged and to an imager platform and provide 
reference points or lines to facilitate the determination of 
the orientation of the patient's head. See, e.g., U.S Patent 

30 No. 4,341,220, which discloses a stereotactic surgical frame 
with fiducial plates that surround the patient's head in the 
fashion of a boxer's headgear and which provides several 
non-collinear fiducial points in cross-sectional scans. Most 
stereotactic frames are fixed to the patient's skull directly, 
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usually by bolts or screws, as noted previously. Clearly such 
methods are not suitable for "routine" office examinations. 

So-called non-invasive, stereotactic devices have also 
been described. The Gill -Thomas stereotactic frame, which is 
5 based on the Brown-Roberts -Wells neurosurgical frame, was 
designed to be used for a series of stereotactic 
radiotherapeutic operations. See, Graham et al . , in 
Radiotherapy and Oncology, (1991) 21:60-62. This device 
requires that it be affixed to the patient by a block, tailored 

10 for individual patients. 

Another device, designed by Laitinen et al., is fixed to 
the patient by means of a nasion support and two ear plugs. 
See, e.g., Laitinen et al., in Surg. Neurol, (1985) 23:559-566 
and U.S. Patent No. 4,617,925. However, this device is then 

15 affixed to the imaging couch or table. Hence, this device is 
able to permit reproducible scans only by relying on the fixed 
position of the patient against the couch or table with respect 
to the machine coordinate system. Again, affixing the patient 
to the machine may makes the patient uncomfortable during the 

20 scan. See, also, U.S. Patent No. 5,330,485, disclosing a 
cerebral instrument guide frame that rests on the bridge of the 
nose (i.e., about the nasion) and which contains plugs for 
insertion into the external ear canals. 

Stereotactic devices are typically fashioned from 

25 precision aluminum alloy and are very expensive for all except 
non-routine use. Moreover, it is usually cumbersome and time 
consuming to affix these devices to the patient, adding to 
their unsuitability for routine examinations. 

The state of the art suggests that stereotactic devices 

30 be equipped with radiographic markers that are visible in scans 
of a patient's head. For example, U.S. Patent No. 4,923,459 
discloses a stereotactic frame that also includes radio-opaque 
rods arranged in the configuration of the letter W N" to 
facilitate localization of a surgical target. U.S. Patent No. 

35 4,608,977 discloses a helmet-like, stereotactic frame that 
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includes such N-shaped "localizing" rod to facilitate the 
determination of the location of a CT scan cross-section. 
Likewise, U.S. Patent No. 4,638,798 discloses a halo-like 
stereotactic frame that has a ring with a plurality of pins of 
5 differing lengths extending therefrom. The relative location 
of a scan can be determined from the relative location of the 
ends of the pins. 

Though such devices can be used to determine the location 
of a head in x-y space, and to determine the relative location 
10 of each imaging "slice," they do not permit the position of a 
head to be fully determined, e.g., as where the head is tilted 
in the imaging plane. 

In addition to the limitations described above, the prior 
techniques are not generally suitable for direct alignment of 
15 images obtained from different imaging modalities. That is, to 
permit the direct comparison of images obtained from different 
imaging modalities, say MR and CT, the patient must be re- 
aligned precisely with respect to the two machine coordinate 
systems. Alternatively, a correction can be made using image 
2 0 processing techniques after a second or subsequent scan has 
been taken. However, image processing has the drawback in that 
the resolution of the processed image is dependent on the 
quality of the scan data set. It would be desirable to alter 
a scan in real time such that scans from different modalities 
25 can be compared directly without the need for image processing. 

In each of the known devices and methods, the anatomical 
coordinates of the patient are fixed in relation to a reference 
coordinate system, that is the machine's coordinate system. 
Thereafter, the machine's coordinate system is used as the 
30 reference coordinate system for each subsequent scan. Because 
of the difficulty in reproducing the machine coordinate system 
or because different machines are invariably associated with 
different, incompatible machine coordinate systems, it has not 
before been possible to relate directly scans from different 
35 imaging modalities. Moreover, it is not always possible to 
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align directly scans from the same imaging modality (e.g., MR 
imagers) when comparing images obtained from machines made by 
different manufacturers. 

It would thus be desirable to have a device and method 
5 whereby the reference coordinate system is independent of the 
machine or imaging modality. It would be desirable, moreover, 
to use a reference coordinate system "personal" to the patient 
as the reference coordinate system and, where possible, have 
the machine's coordinate system fixed or adjusted relative to 

10 that of the patient to provide for scans that are reproducible, 
compatible and superimposable in the same or different imaging 
modalities. A system that enables the taking of imaging scans 
under such a patient reference or "personal" coordinate system 
would be of great utility and would be deemed a significant 

15 advancement in the art. 

In view of the foregoing, it is an object of the invention 
to provide improved devices and methods of non- invasive, 
repetitive, radiographic examination of a subject, particularly 
of the subject's head. 

20 A further object of the invention is to provide such 

devices and methods that are readily amenable for use in 
"routine" examinations, as well as for surgical planning and 
follow-up. 

A still further object of the invention seeks to provide 
25 improved devices and methods of stereotaxis (both invasive and 
non- invasive) , which are low- cost, easy to use, comfortable and 
which provide accurate and reproducible results. 

Yet another object of the invention relates to improving 
methods and apparatuses that can determine fully the position 
3 0 of a head and a scan plane, including when the head is tilted 
in the scan plane. 

Other objects of the invention include providing a way 
or means for comparing directly scans taken by the same or 
different imaging modalities and providing a method for the 
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reproducible placement of external markers, e.g., electrodes, 
on a patient. 

3 . Summary of the Invention 
5 The invention thus provides devices and methods for the 

non-invasive, imaging or radiographic examination of a subject. 
By H imaging* is meant any scanning or spectroscopic technique 
that provides information that can be recorded on a tangible 
medium (e.g., photographic film, slides and the like) or 

10 electronically for storage, later retrieval, or manipulation. 
Moreover, the scanning or spectroscopic technique may also give 
rise to an image that is viewable, e.g., on a screen or 
monitor. The immediate objective, of course, is to provide 
information regarding, or an image of, the internal organs or 

15 tissues of a subject. Such scanning or spectroscopic technique 
or imager can use a wide variety of electromagnetic radiation 
(or for that matter any suitable source of energy) to probe or 
excite internal atoms, ions, molecules, structures, cells, 
tissues, or organs, including but not limited to radio waves, 

20 infrared, ultrasound, ultraviolet, X-rays, electron beam, 
alpha-, beta-, or gamma-rays or particle emissions.. 

Accordingly, the invention provides a stereotactic device 
that is intended for use with an imager and generally comprises 
a frame equipped with localizing means and affixing means. The 

25 localizing means comprises one or more localizing arrays that 
provide one or more imager detectable signals, while the 
affixing means comprises non- invasive fittings for placement 
about the periphery of the subject and which permit the 
reproducible positioning of the frame on the subject . 

3 0 From the signals is derived a personal coordinate system that 
serves as a reference coordinate system for imaging scans taken 
of a subject on which the frame is positioned. This personal 
coordinate system is independent of any machine coordinate 
system. (Contrast the device of Laitinen et al . , for example, 

35 which must be attached to the imaging table or couch to "fix" 
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the device's (and consequently the patient's) coordinate system 
to that of the machine coordinate system.) In a specific 
embodiment of the invention, the localizing means comprises 
localizing arrays or individual reference elements. 
5 Alternatively, the localizing array may be made up of one or 
more reference elements. The resulting three-dimensional 
reference coordinate system is specific or "personal" to the 
subject and is independent of the machine coordinate system. 
The device, as described further below, may be 

10 reproducibly positioned to the subject without the need for an 
invasive affixing means, such as staples, pins and/or bolts 
(i.e., the device of the invention is tt non- invasive" ; non- 
invasive can also mean the absence of a surgical intervention 
or of a breach of a subject's body). 

15 The invention also relates to method of obtaining imaging 

scans of a subject which includes providing a non- invasive 
stereotactic device that is positioned reproducibly on a 
subject and which device establishes a personal coordinate 
system (PCS) associated with the subject. Again, the PCS is 

20 independent of a machine coordinate system (MCS) associated 
with an imager. 

Subsequently, using an imager having an MCS, an imaging 
scan of the subject is taken (including the stereotactic 
device) to establish the PCS of the subject. The MCS of the 

2 5 imager is then manipulated to bring the MCS in substantial 

alignment with the PCS of the subject. One or more additional 
imaging scans of the subject are taken next, with the MCS of 
the imager substantially aligned with the PCS of the subject, 
to provide a first set of imaging scans. 

3 0 In yet another aspect of the invention, a method of 

obtaining imaging scans of a subject taken over different time 
periods is disclosed. The method comprises taking at a first 
time period, using an imager, a baseline imaging scan that is 
relatable to a personal coordinate system (PCS) and a first 
3 5 machine coordinate system (MCS) . The PCS, by definition, can 



WO 97/067< 




PCT/US96/13956 



8 

be regenerated from the subject in a substantially reproducible 
manner independent of the imager's (or machine's) coordinate 
system. At a different time period, using a second imager, at 
least one follow-up imaging scan is taken, which scan is 
5 relatable to the PCS and a second MCS. The second MCS is then 
manipulated, such that the relationship between the second MCS 
and the PCS is substantially the same as the relationship 
between the first MCS and the PCS. At least one additional 
follow-up imaging scan is then taken, which imaging scan can be 

10 superimposed on the baseline imaging scan. The second imager 
may be the same as or different from the initial "baseline" 
imager. It should be understood that follow-up scans can be 
taken over a wide range of time periods, from very short, 
essentially back-to-back scans to much longer time periods of 

15 days, to weeks, to years. 

A further method of the invention relates to yet another 
method of obtaining imaging scans of a subject. The method 
comprises providing a subject with a non-invasive stereotactic 
device that is positioned reproducibly on a subject. The 

20 stereotactic device, when positioned on the subject, 
establishes a personal coordinate system (PCS) associated with 
the subject which, as always, is independent of a machine 
coordinate system (MSC) associated with an imager. The method 
continues with a step of taking, using a first imager having a 

25^ first MCS, at least one imaging scan of the subject including 
the stereotactic device to establish the PCS of the subject and 
to relate the PCS of the subject to the first MCS of the first 
imager. The next step involves taking, using a second imager 
having a second MCS, at least one imaging scan of the subject 

30 including the stereotactic device to reestablish the PCS of the 
subject and to relate the PCS of the subject to the second MCS 
of the second imager. The second MCS is then manipulated, such 
that the PCS is related to the second MCS in substantially the 
same way as the PCS is related to the first MCS. The method 

35 may further comprise taking one or more additional imaging 
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scans of the subject with the second MCS of the second imager 
so manipulated. 

In this manner, the invention provides a device and 
methods by which temporally different imaging scans from the 
5 same or different imaging modalities can be compared directly. 

The device and methods of the invention thereby satisfy 
a long felt need in the art by providing a way to define or 
obtain a personal coordinate system that is independent of the 
machine coordinate system, the time at which the imaging scans 
10 are taken, the operator of the imager, the brand name of the 
imager, the model of the imager, or even the modality of the 
imager. 

An exemplary stereotactic device and method of the present 
invention can thereby facilitate routine examinations of a 

15 patient or subject. In particular, the subject's head can be 
easily and repeatedly examined with the confidence that imaging 
scans taken at different time periods can be superimposed or 
compared directly. In some cases, portions of the patient's 
spine can also be viewed routinely, e.g., during regular, pre- 

20 and post -surgical planning examinations . The imaging modalities 
that can be used to advantage (or used with each other in any 
combination, except that the imaging modality that has a 
"fixed" - not adjustable machine coordinate system should 
preferably be used to obtain the first, baseline, or initial 

25 set of imaging scans) include MR, CT, PET, SPECT, MEG, and 
other such imaging/radiologic scanning or spectroscopic 
techniques . 

These and other aspects of the invention are evident from 
the discussion above and from the more detailed descriptions 
3 0 that follow of the preferred embodiments of the invention. 

4 . Brief Description of the Drawings 

A better understanding of the invention may be attained 
by reference to the drawings, provided herein, in which: 
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Figure 1A depicts a front view of one embodiment of a 
non-invasive, stereotactic device according to the invention; 

Figure IB depicts a top view of one embodiment of a 
non- invasive, stereotactic device according to the invention; 
5 Figure 1C depicts a side view of one embodiment of a 

non- invasive, stereotactic device according to the invention; 

Figure 2 illustrates a technique for determining the exact 
position of reference elements based on a pattern of Mots' 1 
formed in a scan by the device of Figures 1A, IB and 1C; 
10 Figure 3 illustrates another technique for determining 

the position of reference elements based on a pattern of "dots" 
formed in a scan by the device of Figures 1A, IB and 1C; 

Figure 4A depicts a scan of a volunteer's head on which 
is placed the device of Figures 1A, IB and 1C, while Figure 4B 
15 depicts a modified scan rotated and translated according to a 
methodology of the invention. 

Figure 5 is a flow chart showing a diagnostic procedure 
according to the invention for determining the orientation of 
a patient's head. 
20 Figures 6A depicts a scan of volunteer's head on which 

is placed the device of Figures 1A, IB and 1C, while Figure 6B 
depicts another scan of the same volunteer, but by a different 
operator, demonstrating the reproducibility of the inventive 
method. 

25 Figure 7 depicts a volunteer wearing the device of Figures 

1A, IB and 1C. 

Figures 8A-8F depict the specifications for one embodiment 
of the device of the invention. 

30 5 . Detailed Description of the Invention 

The illustrative device and method permit the alignment 
of routine scans to a pre-defined personal coordinate system. 
Once this alignment is performed, the location, orientation and 
size of, e.g., lesions can be reproducibly and precisely 

3 5 determined. It thus becomes possible to compare locations and 
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size from different examinations at different times and sites 
using different machines for a given modality. It is further 
possible using the stereotactic device of the invention to 
provide for imaging scans that can be aligned directly and 
5 compared across different imaging modalities. 

To accomplish further the objectives of the present 
invention, the following detailed description is provided which 
is directed to the lightweight frame, including non-invasive 
affixing means, localizing means and reference elements for 
10 defining a personal coordinate system and a computational 
procedure for the manipulation and/or alignment of coordinate 
systems to permit imaging scans to be taken at different time 
periods which are superimposable or directly comparable. 

15 5.1. Lightweight Frame 



stereotactic device comprises a lightweight frame having 
affixing means comprising non-invasive fittings for placement 
about the periphery of the subject and which permit the 

20 reproducible positioning of the frame on the subject. The non- 
invasive fittings may include, for example, ear fittings and a 
nose fitting. The frame is adapted to fit partially or 
completely around the human head. When present, the nose 
fitting may rest on the nasion or on or about the bridge of the 

25 nose. Also when included, the ear fittings (usually a pair) 
may be inserted into the ears or allowed to rest over them 
(e.g., in substantially the same manner as eye glass frames 
rest on the bridge of the wearer's nose and over the wearer's 
ears) . 

30 Moreover, localizing means, described further below, are 

preferably conveniently positioned in the proximity of the nose 
and ear fittings (e.g., one by the nasion and one for each ear 
fittings for a total of three localizing means) .. In this way, 
a cross-sectional imaging scan using an MR, CT, PET, SPECT, 



In a particular embodiment of the invention / a 
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MEG, and/or such other imager can cut across all three 
localizing arrays of the localizing means in a single slice. 

The frame may be equipped with additional, optional 
features, such as a securing means that facilitates the further 
5 holding or securing of the frame to the subject, especially the 
subject's head. Such securing means include, but are not 
limited to, elastic or inelastic components, e.g., fabrics, 
VELCRO bands and spring assemblies. Preferably, such optional 
securing components are integrated with the aforementioned ear 

10 fittings. Alternatively, these components maybe integral with 
or detachable from the frame. 

Figures 1A - 1C depict front, top and side views of an 
illustrative stereotactic device according to one embodiment 
of the invention. The stereotactic device comprises a curved 

15 frame 10 which is worn by the subject in a manner somewhat 
similar to eyeglass frames. The stereotactic device further 
comprises affixing means, including a nose fitting 12 fixed to 
the midpoint of the frame and a pair of ear fittings 14. In 
some embodiments, the ear fittings are adjustable, slidably 

20 mounted along the frame to accommodate the configuration of the 
subject's head. The nose fitting is preferably placed on the 
nasion, and the ear fittings are preferably placed in, on, or 
over the ears, more preferably in the outermost portion of the 
external auditory canals. 

25 Referring again to Figure 1A, the frame 10 can be of 

variable length, preferably, a length of approximately 15 
inches or so, and is formed into a shape providing for 
comfortable use on a typical human adult. Such a shape may be, 
e.g., a slight curve, a half circle, a "U," or the like. 

30 Different sizes can also be made to accommodate smaller-sized 
heads, such as an adolescent's or child's head. 

Preferably, the materials used for making the frame are 
low-cost, non-magnetic and transparent to the imaging system. 
Moreover, the materials should be durable and preferably 

35 amenable to repeated sterilization (especially when 
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implementing non-disposable frames) . Preferred materials for 
the frames include relatively rigid thermoplastic materials, 
including most synthetic polymers but most preferably 
plexiglass . 

5 

5.2. Localizing Means 

In this aspect of the device, the frame is equipped 
with localizing means comprising one or more, preferably at 
least three, of what are referred to herein as localizing 
10 arrays. Each array provides one or more imager detectable 
signals from which the personal coordinate system of the 
subject is derived. With the personal coordinate system 
serving as the reference coordinate system for the initial, 
baseline imaging scans and each subsequent follow-up imaging 
15 scans, the stereotactic examination of the subject is 
facilitated . 

In a particular embodiment of the invention, the 
localizing array is made up of one or more reference elements, 
preferably including a pair of reference elements. A suitable 

20 reference element may comprise an elongate component, e.g., a 
cylindrical rod or tube. Other localizing means and 
corresponding localizing arrays would be apparent to one of 
ordinary skill on appreciation of the disclosure provided 
herein. For example, a localizing array may comprise a 

25 supporting means, such as a cylindrical guide directed toward 
the body of the subject, into which a digitizing w pen" can be 
inserted. The digitizer can, in turn, emit signals that can be 
observed and/or recorded on an imaging scan. Also, depending 
on the number of localizing arrays, which are provided with a 

30 frame, the reference elements can include spherically shaped, 
egg-shaped, or irregularly shaped "opaque" (more, below) 
materials to create points in space. 

The preferred elongated pair of reference elements can 
have a number of configurations but are preferably arranged in 
5 a spaced-apart n X" configuration. In one embodiment of the 
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invention, the elongate components of each pair form an tt X" but 
do not intersect. Rather, the elongate components (e.g., rods 
or tubes) are staggered from one another so as to be spaced 
apart. Thus, e.g., the reference elements of an array, 
5 positioned in the proximity of, or at, the nasion, inscribe an 
W X" configuration when viewed initially from the front of the 
head and form spaced apart, substantially parallel lines when 
viewed from the top of the head at approximately a ninety 
degree angle from the initial view. Generally, the reference 

10 elements can *be spaced apart by any practical distance, but 
preferably range from about 0.1 to about 1 inch, more 
preferably, about 0.2 to about 0.5 inch. 

The reference elements are constructed or filled with 
materials that produce distinctive features in an imaging scan, 

15 such as a radiographic scan using an MR imager or CT scanner. 
For example, the reference elements can be filled with doped 
water, which is relatively "opaque" (i.e., give rise to 
detectable signals in the course of the imaging scan) to MR, or 
they can be constructed of a material that is itself opaque to 

20 MR (e.g., stainless steel) . When such a stereotactic device is 
worn on, e.g., the head of, a patient being scanned, the 
reference elements appear as Mots" (in the case of 
cylindrical, rod* like, or spherical components) or other 
distinctive features (in the case of non-cylindrical or 

2 5 irregularly shaped components) in the resulting scan. 

Referring again to Figures 1A, IB and 1C, a localizing 
means is shown comprising three localizing arrays, each array 
in turn comprising paired reference elements 16A/16B, 18A/18B 
and 20A/20B. The three localizing arrays are each mounted on 

30 three ear and nose fittings, as shown. Each pair of reference 
elements is arranged in the form of an *X," when viewed from 
the front (tubes 16A/16B) or side (tubes 18A/18B and 20A/20B) 
and is constructed from, or preferably filled with, a 
radiographically detectable substance, such as doped water. 
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In cross-sectional, radiological scans, such as MR images, 
the paired reference element tubes appear as two dots at each 
fitting or location, six dots in toto. Those dots are shown in 
Figures 4A and 4B, where dots a and b correspond to the cross - 
5 sectional image of the paired reference elements 16A/16B, 
respectively; dots c and d correspond to the image of the 
paired reference elements 20A/20B; and dots e and f correspond 
to the image of the paired reference elements 18A/18B. 



10 a variety of materials or combinations thereof, including 
plexiglass or other non-magnetic, radio -transparent material 
that is filled with a radio-opaque substance. Alternatively, 
these elements may be formed from a radio-opaque substance, 
such as steel or other metal. The reference elements can be 

15 circular in cross-section (e.g., with a radius of approximately 
0.05 to approximately 0.4 inch, preferably about 0.2 inch, more 
preferably about 0.1 inch) or of any other cross-sectional 
shape readily discernible on a scan. Where the reference 
elements are in the form of cylindrical tubes, they are 

20 generally from about one to about five inches in length, 
preferably about two inches in length. 

In one embodiment of the invention, in which the reference 
elements are paired, the members of each pair are generally 
spaced apart from one another. This spacing may be about 0.1 

25 to about 0.5 inch but is most preferably about 0.2 inch. 

The reference elements are mounted on plexiglass or other 
non-magnetic, radio- transparent support of any size suitable 
for providing a supporting surface for the tubes without making 
the stereotactic device unwieldy. For example, squares of 

30 approximately one-by-one inch to five-by-five inches, 
preferably about two-by-two inches, more preferably about 1.4 
x 1.4 inches, are suitable for use as the support for the 
reference elements. Hence, in a particular embodiment of the 
invention, the combination of support and the reference 

35 elements comprise a localizing array. In another embodiment of 



As stated above, the reference elements may be made of 
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the invention, a support may also serve as a guide for a 
digitizing "pen" that can be manipulated by a physician or 
technician and which generates detectable signals at locations 
that can be observed and pinpointed in an imaging scan. In 
5 still other embodiments of the invention, supports may be 
dispensed with and the reference elements are found directly on 
or in the frame (e.g., the elements may be an integral part of 
the frame) . 

As used herein the term "radio-opaque" and 
10 "radiographically opaque" refer to materials, such as doped 
water, which are visible on an MR I scan or on such other 
radiologic scan. Of course, a generally "opaque" material can 
be chosen so that detectable signals can be observed in each, 
or some, of MR, PET, SPECT, CT, MEG, or X-ray, as the case may 
15 be, and using which material portions of the stereotactic 
device (e.g., the reference elements) according to the 
invention are made. Likewise, the term "radio- transparent" 
refers to materials, such as plexiglass, plastics and most 
synthetic polymer materials, which are generally not visible on 
an MR I scan or on such other radiologic scan (e.g., PET or 
SPECT), with which portions of the stereotactic device (e.g., 
the frame) according to the invention is used. Such radio 
transparent materials can also be used in combination with 
radio opaque materials, as would be apparent to one of ordinary 
skill. It should be noted that plexiglass, plastics and most 
synthetic polymers are visible, and thus opaque, in certain 
imaging modalities, such as CT or X-ray. 

5.3. Imaging Procedure 

To analyze a scan that includes an image of the 
reference element tubes, a reference point is first defined; 
most suitably, the cross-point of each pair of reference 
element tubes in the above -described n X" configuration can 
serve as this reference point. Assuming the image is in the 
in-scan plane, e.g. as in Figures 4A and 4B, the x and y 
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coordinates of the reference point can be found as the mid- 
point between two dots on the image. The z coordinate, or the 
through-plane position, is determined from the distance of the 
two Mots" generated by each tube. As a result, the three- 
5 dimensional positions of three reference points (one from each 
pair of tubes) with respect to the MCS are obtained from a 
single axial image . These reference points, in turn, define 
the personal coordinate system. 

The machine coordinates (i.e., the machine coordinate 

10 system) can then be adjusted to rotate and translate the scan 
to a predetermined plane defined by the personal coordinate 
system. In this manner, the machine coordinates system can 
always be made to coincide with reference points that are 
independent of the machine and which are determined in relation 

15 to the specific patient or subject. Hence, a series of imaging 
scans can be obtained of the same patient at different times, 
independent of the operator or the specific imaging instrument, 
all of which are substantially superimposable on one another, 
including the initial set or baseline set of imaging scans. 

20 Also, images from differing imaging modalities can be "merged" 
by using the stereotactic devices and methods of the invention 
to provide a composite image comprising superimposed images 
taken from one or more different imaging modalities. 

More particularly, analysis of an imaging scan to 

25 determine the position of the reference elements of the 
aforementioned stereotactic device can be accomplished by the 
methods of the invention, such as those described in detail, 
below. 

5.3.1. Two- Scan Exact Positioning Method 
3 0 From two or more imaging slices that cut 

through both tubes of each of the three localizing arrays, the 
exact location of the reference point can be calculated 
mathematically as follows. 

Referring again to Figure 2, points a, b, c and d are 
3 5 where the two planes of the two imaging slices cut through the 



WO 97/067- 




PCT/US96/13956 



18 

staggered n X" (thin lines) . The x, y and z coordinates in the 
initial imaging coordinate system can be read from the images. 
In addition, e and f are points on the tubes where the two 
tubes are closest to each other, and the reference point, g, 
5 can be defined as the mid-point between e and f . The objective 
is to find the coordinates of g, given coordinates of a, b, c 
and d. 

Two unit vectors, parallel to each of the tubes, V x and V 3 , 
can be defined as: 

10 

b-a 



15 

d-c 

20 

Also, a third vector which is orthogonal to vl and v3 can be 
defined as: 



V 2 = V, x V 3 

25 

where x denotes the vector (outer or cross) product. 
Then line segments a-e-f-d can be expressed as: 



AV X + BV 2 + CV 3 

where A, B and C are (scalar) distances between a and e, e and 
f, and f and d, respectively. 

Since the vector a->e + e->f + f->d is also equal to a->d, 
then V 4 can be defined as: 



V 4 = d - a 
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and 



v 4 - tv^v,] 



Since V x , V 2 and V 3 are orthogonal, A, B and C can be 
obtained as: 



where "T" is transpose. 



20 



25 



30 



35 



The point g is: 
g s a + AV X 



B 
2 



_ v 2 



For the above-described computational method to work well, 
the relative position of the two reference element tubes in 
each localizing array must be defined by some known, non-zero 
angle. The above steps are simplified by assuming that the two 
tubes in each localizing array are orthogonal to each other. 
However, the above -described method also works as long as two 
tubes have any pre-determined, non-zero angle. In that case, 
however, the inversion of [V^Vj] becomes necessary, as will be 
apparent to one skilled in the art. 

It is further preferable for the sake of simplicity in 
calculating the scan position to have the angles between the 
tubes in each localizing array to be the same. Again, however, 
this is not critical under the practice of the invention so 
long as the angles of the tubes in each localizing array are 
known . 



40 
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5.3.2. Single-Scan Positioning Method 

It is also possible to determine the 
position of the reference point from a single slice relatively 
accurately, although mathematically not exact. Assuming the 
5 imaging scan is in the x-y plane, the x and y coordinates of 
the reference point can be found as the mid-point between two 
dots on the image. The z coordinate, or the through-plane 
position, can be determined from the distance between the two 
dots. Since the two tubes are arranged in a staggered U X" 
10 configuration, the direction in the z axis can be determined 
knowing the angle between the elongated components. As a 
result, the positions of three reference points in all three 
dimensions are obtained from a single social image. 

Referring now to Figure 3 , points a and b are where the 
15 imaging plane (thin line) cut through the staggered W X" (thick 
lines) . The distance D is defined as the distance between the 
points a and b on the image and corresponds to the distance 
between the two tubes that make up the W X" of the localizing 
array. The x and y coordinates of the reference point, m, can 
20 be found on the image as the midpoint between the dots. The z 
coordinate of the reference point, z, can be calculated as: 

•"/(SI'-U)' 

5.3.3. Computational Method 
The detailed calculation method for 

determining the orientation of a patient's head is as follows: 
Let pi, p2 and p3 be vectors representing three reference 
points in the initial image coordinate system corresponding to 
a/b, c/d and e/f in Fig. 4A, respectively. These points are 
calculated as described above. Then, unit vectors representing 
the new, desired personal coordinate system, x, y, z, can be 
calculated as: 
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10 



15 
20 



25 



30 



35 



x = p3 - p2 

z = x x (pi - p2) 

y = x x z 

where x denotes the vector product operation. A rotation matrix 
that rotates these into xyz axes is obtained as: 



R 




where n T n is transpose and the translation vector t, is 



A general diagnostic method according to the invention is 
outlined in Figure 5, including determining a standard 
reference coordinate system that is based on the personal 
coordinate system. First, a stereotactic device, as described 
above, is placed on the patient. Next, the patient is placed 
in the imaging device while wearing the stereotactic device. 
The patient position is adjusted so that an axial image cuts 
across all tubes in at least three localizing arrays. A scan 
is taken subsequently. An exemplary resulting scan is shown in 
Figure 4A. 

From such a scan, the three patient reference points of 
the personal coordinate system are calculated as described 
above. From these points, a rotation matrix and a translation 
vector are calculated to adjust the initial imaging coordinate 
system to the personal coordinate system, as defined by the 
stereotactic device. On the MR imager, these transformations 
are applied to the gradient fields, receiver frequency and 
phase to change the machine coordinate system. 

Next, another axial imaging scan is obtained in the new 
standard reference coordinate system and the Mots" are checked 
to ensure that the Mots" are aligned at the predetermined 
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positions. The scan of Figure 4A # modified in accord with this 
method of the invention, is shown in Figure 4B. 

Once a standard reference coordinate system is determined 
by the foregoing steps, then a full set of examinations is 
5 obtained using the defined rotation matrix and the translation 
vector. The entire image alignment process can be performed in 
a minimal time as compared to previously known scanning 
techniques in which an attempt is made to adjust or fix the 
patient coordinate system each time relative to the coordinate 
10 system of the imager. 

5.4 . Additional Supporting Disclosure 

A still further aspect of the invention provides a 
stereotactic device for non-invasive stereotactic examination, 

15 particularly of the head of a subject, comprising a frame that 
is reproducibly positioned on the subject. The frame may have 
any number of localizing means comprising radio-opaque 
reference elements to provide for a multi -dimensional reference 
coordinate system. In a preferred embodiment, six reference 

20 elements are included, which are arranged in pairs to provide 
three localizing arrays. 

The frame of the stereotactic device is preferably made 
of a non- radiographic material, preferably plexiglass. Other 
suitable radio- transparent materials include, but are not 

25 limited to plastics, synthetic polymers, or other carbon-based 
materials of some structural rigidity, such as poster board, 
cardboard, or even graphite. 

In a particular embodiment, the frame provides a three- 
dimensional framework for the at least four reference elements 

30 to provide a three-dimensional personal coordinate system. In 
such an embodiment, the localizing means comprises four or more 
localizing arrays each comprising a reference element. Each 
reference element, in turn, defines a point in space, three of 
which points define a unique plane and the fourth point lying 

35 outside the unique plane. 
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In yet another embodiment three reference elements can be 
designed into the device to provide inherently the fourth 
reference element necessary to define a three-dimensional space 
(e.g., as described earlier, confining two of the three 
5 reference elements to a pre-determined plane) . 

The reference elements are positioned on the frame at 
predetermined positions to provide the three dimensional 
coordinate system. As described above, the reference elements 
can be paired in an orthogonal arrangement and be constructed 
10 or filled with materials that produce Mots" or other 
distinctive features in an imaging scan plane. These elements 
make up the localizing means that provides a reproducible 
localization of the multi-dimensional personal coordinate 
system. Thus, the reference elements may include markings such 
15 as small paint spots, indentations, slots, grooves, and the 
like, in or on the frame whose positions can be entered into or 
recorded in a scan via, e.g., an MRI- compatible digitizer. 

Hence, in a preferred embodiment of the invention a 
stereotactic device for use with an imager comprises a frame 

2 0 equipped with localizing means and affixing means. The 

localizing means comprises three or more localizing arrays each 
equipped with two or more reference elements that together 
provide six or more imager detectable signals. The affixing 
means comprising three or more non- invasive fittings for 

25 placement on the subject's nasion and in or about the subject's 
ears. The fittings, together, permit the reproducible 
positioning of the frame on the subject and the reference 
elements each provide an imager detectable signal .' The personal 
coordinate system is then derived from these signals, 

30 collectively.. 

The reference elements preferably comprise elongate 
components. For example, two elongate components are arranged 
in an n X" configuration and are related to one another by a 
predetermined angle, e.g., about 90 degrees. Further, the 

3 5 reference elements comprise radiographically opaque or semi- 



WO 97/067< 




PCT/US96/13956 



24 

opaque material, including steel or doped water. The signals 
obtained from the reference elements are preferably detectable 
by MR, CT, PET, SPECT, EEG, or MEG. In an alternative 
embodiment, the imager detectable signals are provided by a 
5 digitizer used in conjunction with the one or more localizing 
arrays. 

The preferred device may further comprise a securing means 
that facilitates the securing of the frame to the subject. 
Such securing means may include an inelastic or elastic 
10 component. 

In use, the stereotactic device is reproducibly positioned 
on the head of a subject. The position of the reference 
elements with respect to the desired location of the subject is 
then digitally recorded, e.g., with an MRI -compatible 

15 digitizer. The digitized positions of the reference elements 
provides a reproducible personal coordinate system. 
Subsequently, the positional rotation and translation 
necessary to bring the machine coordinate system into alignment 
with the personal coordinate system can be determined, 

20 preferably by a digital processor. After the machine 
coordinate system is adjusted to align with the personal 
coordinate system, one or more radiographic scans are taken. 

In subsequent scans, the patient again reproducibly dons 
the stereotactic device and the positions of the reference 

25 elements are digitally recorded. The second imager is then 
adjusted to the personal coordinate system by means of a 
translation and/or rotation of the second imager's machine 
coordinate system, as above. The second and subsequent scans 
are then taken, again in the patient coordinate system. By 

30 scanning each time in the patient coordinate reference system, 
all scans can be directly compared. 

Described above are improved devices and procedures for 
non-invasive, radiographic analysis, particularly stereotactic 
head examinations, e.g. , in connection with stereotaxy or other 

35 similar surgical procedures, meeting the objects of the 
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invention. Some of the more noteworthy features of the 
invention include, but are not limited to: 

1. A stereotactic device that can be made with extremely 
low cost and which can be made in disposable or re-usable form; 
5 2. A method of aligning the imaging plane in which only 

one axial image through the device is necessary to perform the 
subsequent alignment; the only human intervention required is 
the identification of the reference points (Mots") of the 
personal coordinate system on the localizing axial image; the 
10 rest of the process can be implemented in the imagery and 
processed automatically; moreover, with an automated mechanism 
to identify the "dots," the entire process can be fully 
automated; 

3 . Imaging scans acquired according to the invention are 
15 alignable to a reproducible reference coordinate system, so 

that the images from different examinations can be compared 
directly; these examinations could be performed at different 
institutions, with different imagers, etc., as long as all 
examinations are taken with a device of the invention 
20 positioned on the subject; 

4. A device and method best used in obtaining MR images; 
however, other imaging modalities, such as CT, PET, SPECT, or 
MEG, are also applicable; these other imaging modalities can 
then be cross-referenced with MR imaging scans. 

25 The stereotactic device of the invention is exemplified 

as an eyeglass-like structure. However, the invention may be 
of any suitable structure that can support the localizing 
arrays or reference elements and which can be reproducibly 
positioned on the subject. Hence, the invention provides a way 

30 of mapping an anatomical region of a subject which can be 
related to a personal coordinate system that is independent of 
the machine coordinate system. Furthermore, the device and 
methods of the invention may also be amenable to veterinary 
applications . 
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It should thus be apparent that the present invention 
provides a method of obtaining imaging scans of a subject 
comprising the steps of: providing a subject with a non- 
invasive stereotactic device that is positioned reproducibly on 
5 the subject and which establishes a personal coordinate system 
(PCS) . The PCS is associated with the subject independent of 
a machine coordinate system (MSC) that is associated with an 
imager; taking, using an imager having an MCS, an imaging scan 
of the subject including the stereotactic device to establish 

10 the PCS of the subject; manipulating the MCS of the imager to 
bring the MCS in substantial alignment with the PCS of the 
subject; and taking one or more additional imaging scans of the 
subject with the MCS of the imager substantially aligned with 
the PCS of the subject, to obtain a first set of imaging scans. 

15 The method of the invention may further comprise repeating 

the above-mentioned steps at a second time period, using a 
second imager, to obtain a second set of imaging scans. 
Subsequently, at least one imaging scan of the first set can be 
compared with at least one imaging scan of the second set. In 

20 this way, an operator has the opportunity to note and make a 
record of any previously undetected anatomical feature of the 
subject. Moreover, observations can be made of any changes in 
any previously detected anatomical feature of the subject. As 
mentioned above, such anatomical features may be of anything 

25 that can be of interest to the subject or the medical 
practitioner, including but not limited to lesions, tumors, or 
features that may indicate a pathological condition. 
Adventitiously, the stereotactic device is positioned 
reproducibly on the subject's head. 

30 The second time period of the disclosed method represents 

an elapsed time from the taking of the first set of imaging 
scans to the taking of the second set of imaging scans. This 
elapsed time may, of course, be any time period appropriate to 
the examination process, including periods that are very short 

35 (e.g., essentially back-to-back scans), intermediate (e.g., 
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days to weeks) , or very long (e.g., years) . For instance, the 
elapsed time period may be about 15 to about 45 minutes, a few 
hours, one day to about one week, about one week to about one 
month, about one month to about six months and about six months 
5 to about one year. This second time period may even represent 
an elapsed time from the taking of the first set of imaging 
scans to the taking of the second set of imaging scans of about 
one year to about five years. 

Separately, the invention provides a method of obtaining 

10 imaging scans of a subject comprising the steps of (a) 
providing a subject with a non- invasive stereotactic device 
that is positioned reproducibly on a subject and which 
establishes a personal coordinate system (PCS) associated with 
the subject which is independent of a machine coordinate system 

15 (MSC) associated with an imager; (b) taking, using a first 
imager having a first MCS, at least one imaging scan of the 
subject including the stereotactic device to establish the PCS 
of the subject and to relate the PCS of the subject to the 
first MCS of the first imager; (c) taking, using a second 

20 imager having a second MCS, at least one imaging scan of the 
subject including the stereotactic device to reestablish the 
PCS of the subject and to relate the PCS of the subject to the 
second MCS of the second imager; and (d) manipulating the 
second MCS, such that the PCS is related to the second MCS in 

25 substantially the same way as the PCS is related to the first 
MCS. 

In this method of the invention, the second MCS is 
substantially aligned with the first MCS. In a specific 
embodiment, the first imager may be an imager in which the 

3 0 machine coordinate system cannot be adjusted to the personal 
coordinate system (i.e., the machine coordinate system is 
fixed, as with a CT scanner) . If so, the second imager is one 
whose machine coordinate system is adjustable, preferably, an 
MR imager. The method involving different imaging modalities 

3 5 may further comprise forming a composite image including 



WO 97/067< 




PCT/US96/13956 



28 

information from at least one imaging scan taken using the 
first imager and information from at least one imaging scan 
taken using the second imager. 

In addition, the present invention also relates to a 
5 method of spatially aligning at least two radiographic imaging 
scans of the head of a subject taken by a radiographic scanning 
device comprising: (a) radiographically scanning the head of a 
subject wearing a stereotactic device of the invention to 
provide a first radiographic imaging scan containing a 

10 plurality of* reference points defining a first personal 
coordinate system; (b) relating the first personal coordinate 
system with a personal coordinate system obtained from a prior 
radiographic imaging scan of the subject reproducibly wearing 
the stereotactic device; (c) adjusting the radiographic 

15 scanning device to align the first and prior personal 
coordinate systems; and (d) obtaining one or more additional 
radiographic imaging scans of the subject with the radiographic 
scanning device so adjusted. 

The following additional examples are provided to further 

20 illustrate preferred aspects of the invention. Nothing in 
these additional examples should be construed to limit the 
invention in any way. 

6 . Examples 

25 

6.1. General Method 

In a generic embodiment of the invention, a pulse 
sequence to rotate the imaging coordinate system is executed on 
a clinical MR imager (e.g., a 1.5T Signa, GE Medical Systems) . 
30 The rotation matrix is determined by a conventional digital 
data processor (e.g., a personal computer) programmed to 
calculate the rotation matrix from the positions of the 
reference points of the frame. A listing of the relevant 
computer programs are attached hereto, as Appendices A and B. 
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Those skilled in the art will appreciate that rotation of 
the imaging coordinate system and determination of the rotation 
matrix can be implemented on the image as a single step. In 
that case, the necessary operation is reduced to identifying 
5 the positions of the reference points (e.g., 6 points in a 
specific embodiment) on the first image. Hence, the required 
sophistication of the technologist is greatly reduced by the 
inventive method. 

10 6.2. Use* of an MRI -Compatible Digitizer 

In a specific embodiment of the invention, the 
reference points in the first scan are entered into a digital 
data processor by means of an MRI -compatible digitizer. In 
this embodiment, the technologist need only mark the position 

15 of the reference elements on the first scan and the digital 
data processor calculates the patient-defined personal 
coordinate system from the digitized entries using, e.g., the 
computational examples given above. 

20 6.3. Use of Electrodes for Determining Brain 
Structures 



In yet a further embodiment of the invention, the 
25 correlation of electrical activity in the brain with specific 
brain structures may be reproducibly accomplished. In this 
embodiment, an MR imaging scan is taken of the patient wearing 
the stereotactic device to obtain a "baseline" scan. (It 
should be noted that as long as the device of the invention is 
30 worn, the baseline scan may be taken before or after the EEG.) 

An EEG is then taken of the patient, with the device in 
place. One or more electrodes, typically 15-30 electrodes, are 
. then positioned on the head of the patient. The position of 
each of the electrodes is then entered into a digital signal 
35 processor, e.g., as described above, to define the electrode 
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positions with respect to the patient coordinate system as 
defined by the device. 

If any aberrant, unusual, or abnormal traces are observed 
from the EEG, the identity and position of the electrode or 
5 electrodes giving rise to such traces can be determined. The 
MR imager is subsequently utilized to provide an imaging scan 
across the positions of the electrodes of interest. One or 
more scans are taken to obtain images of the patient's anatomy, 
which may give rise to the aberrant traces. Through this 
10 alignment technique, electrical signals detected by the 
electrodes can be correlated with particular brain structures 
to a high degree of accuracy and reproducibility. 

6.4. Follow-Up Electroencephalography (EEG) 

15 The general positions of each electrode can be 

reproduced, if desired, to correspond with an initial EEG, 
using the device of the invention. Alternatively, each new 
configuration of the electrodes can be recorded digitally 
relative to the PCS established by the device, and MR images 

20 can be taken selectively based thereon. 

When seeking a reproducible placement of electrodes for 
use in electroencephalography (EEG) , a series of MR imaging 
scans is taken of a patient wearing the stereotactic device of 
the invention, as described above, to establish the personal 

25 coordinate system. One or more electrodes are then placed on 
the patient's head, and the position of the electrodes is 
digitized, as described above. The location of the electrodes 
is thereby determined and recorded relative to the personal 
coordinate system . 

30 Upon subsequent EEG examinations, the patient wears the 

stereotactic device of the invention. One or more scans are 
taken to reestablish the personal coordinate system. One or 
more electrodes are again placed on the patient so that their 
positions correspond to their original positions in the initial 

35 EEG/MRI scan. Correct placement of the electrodes is checked 
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by comparing the prior scan and electrode placement with that 
of the subsequent scan. 

Accordingly, any changes in the electrical activity of the 
patient can be attributed to specific electrodes and/or patient 
5 anatomical features as revealed by the MRI . 

6.5- Imaging Across Different Modalities 

A further aspect of the invention provides for the 
direct comparison of imaging scans taken with different 

10 modalities, particularly wherein at least one of the modalities 
has a coordinate system that is not adjustable. In this 
embodiment, e.g., a first scan having fixed machine coordinates 
is taken of the subject wearing the stereotactic device of the 
invention. The reference elements are located in the scan and 

15 a personal coordinate system is established therefrom. 
Additionally calculated is the translation and rotation of the 
machine coordinate system required to align the machine with 
the personal coordinate system. 

A second scan is later taken with a modality in which the 

20 machine coordinate system can be adjusted (e.g., MR imager) 
with the patient reproducibly wearing the stereotactic device. 
The reference elements are identified in the second scan and 
the personal coordinate system determined. The value of the 
difference in the personal coordinate system and the machine 

25 coordinate systems for the initial and second scans are 
determined. 

With the common personal coordinate system in place, the 
second imager is then adjusted to translate and/or rotate the 
second imager coordinate system so that the difference between 
30 the second image coordinate system and the second personal 
coordinate system coincides with the difference between the 
first image coordinate system and the first personal coordinate 
system. 

Subsequent scans are then taken with the second imager 
35 aligned with the scan plane of the first imager, permitting a 
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highly reproducible direct alignment of scans taken in the two 
different modalities based upon the common personal coordinate 
system. 

The stereotactic device and procedure of the invention can 
5 thus be used to align images from modalities such as computed 
tomography (CT) , positron emission tomography (PET) , single 
photon emission computed tomography (SPECT) , electro- 
encephalography (EEG) or magnetoencephalography (MEG) , and the 
like. 

10 In the case of CT and the device of the invention, an 

illustrative procedure is as follows (NB: It is not possible to 
freely position the CT slice, hence, MR images have to be 
rotated to the CT plane.): 

First, the CT examination is done with the device on. The 

15 position of the reference points (mid-point of the X-shaped 
rods) in the CT coordinate system is calculated. Next, an MR 
axial slice through the device is taken, and the MR coordinate 
system is rotated/translated to coincide with the CT coordinate 
system, i.e., the second machine coordinate system is 

20 substantially aligned with the first machine coordinate system. 
Imaging slices that are of the same location/thickness as those 
of the CT examination are then taken. 

As is readily evident to those skilled in the art, similar 
procedures can be employed for aligning other imaging 

25 modalities. In this way, all imaging capabilities can be 
linked through the use of MR under the practice of the 
invention. 

6.6. Process for Determining Subject Head Orientation 
30 Still another aspect of the invention provides a 

method for determining the orientation of a subject's head. 
The procedure includes placing a plurality of, preferably 
three, localizing arrays of the "X"-type, described above, 
about the head of a subject patient, e.g., substantially 
3 5 adjacent the nasion and ears, or in other locations spaced 
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apart about the circumference of the head. The procedure 
further includes taking a cross-sectional imaging scan of the 
head in the vicinity of the localizing arrays and determining 
the orientation of the head based on the pattern of Mots" (or 
other distinctive features generated by the reference elements) 
in the imaging scan. This operation may be accomplished 
conveniently through use of the computer- implemented program 
appended hereto. Adjustments to the orientation of the 
subject's head can then be performed as needed or desired. 



6.7. Further Examples 

The stereotactic device and the correction procedure, 
described above, are used with both phantoms and human 
volunteers. The standard quadrature head coil is used for 
15 taking the imaging scans. 

Figures 4A-4B show a typical set of imaging scans, which 
is acquired with a volunteer. As noted above, Figure 4A is the 
initial axial image taken through the device, and Figure 4B is 
the image taken after the manipulation step, e.g. , the rotation 
20 and translation steps, is applied. All three pairs of dots are 
aligned with each other, indicating the imaging plane is now 
cutting through the three reference points. Also, the 
reference points are at the pre-determined locations. 

To demonstrate the reproducibility of the method, the 
25 following study is carried out: First, the device is positioned 
on the volunteer, and the first axial MR slice is taken. 
Rotation/translation correction is applied and a set of 
multi-slice images is taken. The volunteer is withdrawn from 
the magnet. The device is removed, and replaced again on the 
30 subject by a different technologist. Images are taken again in 
the same manner as in the first scan set. 

Several slices from these two examinations are shown in 
Figures 6A and 6B. Figure 6A shows the images taken by the 
first operator, and Figure 6B shows the images taken by the 
35 second operator. (Imaging parameters are: modified RARE 
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sequence, 30 cm FOV, 256 x 256, 3 mm slice, single slice axial 
scan (Figures 4A and 4B) and multi-slice axial scan with 20 cm 
FOV (Figures 6A and 6B) . ) 

As is evident from these pictures, a high degree of 
5 reproducibility in scan sections is provided by the device and 
method of the invention. 

The accuracy in reproducing slice position and orientation 
between different time points is evaluated in phantom and human 
volunteer studies. The error in phantom studies, when using 

10 the mathematically exact, two- scan method, is about 1 mm or 
less! In human volunteer studies, the error is somewhat higher 
due to patient motion, though still readily acceptable for use 
in comparison of the sizes and locations of lesions from one 
exam to another. In particular, the device of the invention 

15 allows for routine scanning, in part because of the little time 
(2 to 3 minutes) required for device positioning, reference 
scanning and determination of reference points. 

Those skilled in the art 'will appreciate that the 
embodiments described above are exemplary and that other 

20 embodiments incorporating alterations and modifications therein 
fall within the scope and spirit of the invention. Each of the 
references mentioned above is incorporated by reference herein. 
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Appendix A 

° Brigham & Womens Hospital 1995 

calculate rotation matrix for the "frame" 
find the rotate matrix and the translate vector 

2-21-1995 K. Oshio 

last modified on 8-1-1995 (correction step) 

(localizer : 1 axial plane, target : axial plane) 

<- 0 was : -> 1 

-> 1 was : <- 2 



2v A 3 4 A v5 
— plan — 

spinor output (for typing convenience) 

•/ 

#include <stdio.h> 
#include <math.h> 

main() 
{ 

double x[6]. y[6]; /• ref points in ref plane •/ 

double p0[3], pi (3], p2[3], p3[3]; /* ref vectors •/ 

double X[3], Y[3]. Z[3]; /• axes •/ 

double R[3][3]; /• rot matrix •/ 

int i; 

char corr[256]; 

double thap, thlr; 

/* (1) input 6 points •/ 

get _pointfref 0 (x y) ? \ &x[0], &y[0]); 
get_pointCref 1 (x y) ? &x[l), &y[l]); 
get joimTref 2 (x y) ? \ &x[2], &y[2]); 
get joimTref 3 (x y) ? " f ftx[3], &y[3]); 
get_point( H ref 4 (x y) ? \ &x[4], &y[4]); 
getj>oint("ref 5 (x y) ? &x[5], &y[5]); 

/• (2) convert to 3 ref vectors •/ 
/• x */ 

pl[0]«(x[0] + x[l])/2; 
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p2[0] = <x[2] + x[3]) / 2: 
p3[0) = (x[4] + x[5]) / 2; 
/• y •/ 

pi{i] = (y(0] + y[n)/2; 

5 P 2{1] = (y[2] + y[3]) / 2: 

p3[l] = (y[4] + y[5]) / 2; 

/• 2 •/ 

distjo 2(&pl[2], x[0], y[0], x[l], y[l]); 

dist_to_z(&p2[2), y[2], x[2], y[3], x[3]); 

10 dist_to_z(&p3[2], y(5], x[5], y(4], x[4]); 

print_vectorC*pr, pi); 
print_vector("p2",- p2); 
print_vector("p3". p3): 

15 

/• (2.5) correction "/ 

printfC'correciion step ? (y/n) "); 
scanf("%s". corT); 
if(corr[0] — V) { 
20 thap = atan(((p2[2] - p3[2]) / 2 - p I (2)) 

/«p2[l] + p3[l])/2-pl[l])); 
thlr = atan((p3[2] - p2(2]) / (p3[0] - p2{0])); 
pl[0] -= pip] • sin(thlr); 
pi (2] '= cos(thlr) * cos(thap); 
25 p2[l] += p2[2] * sin(thap); 

p2[2] '= cos(thlr) • cos(thap); 
p3(l] +=p3[2] • sin(thap): 
p3[2] •= cos(ihlr) • cos(thap); 

30 printf("thap = %f thlr = %f\n\ thap, thlr); 

print_veaor("pr, pi); 
print_vector("p2". p2); 
print vectorCpS", p3); 

} 

35 

/• (3) make 3 axes •/ 

vsub(X. p3, p2); /* X V 

vsub(Y, pi, p2); I* temp V 

vemul(Z, X, Y); I* Z •/ 

40 vemul(Y. Z. X); /* Y •/ 

vnorm(X); vnorm(Y); vnorm(Z); 

/• (4) rotate mat V 

for (i - 0; i < 3; H-+) { 
45 R[0][i] = X[i]; 

R[l][i] - Y[i]; 
R(2][i] = Z[iJ; 
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} 

vmmuipc R, X); 
vmmui(Y, R. Y); 
vmmul(Z, R. Z); 
vadd(p0, p2, p3); 
vcmul(p0, 0.5, pO); 

/* print result •/ 

mtrans(R); /* inverse •/ 

vctnul(pO, -1.0, pO); /* inverse V 

/• correction •/ 

p0[2] - .p0[2]; /• I+. S- V 

R[0][1] = -R[0][l]; 
■R[1][0] - -RIipj; 
R[0][2] = .R[0][2]; 
R[2][0] = -R[2][0J; 

printf( M \n M ); 

prim_matrixCR : \ R); 
print_vectorCp0'\ pO); 



dist_to_z(z. xh yl, x2, y2) 
double *z; 

double xl, x2, yl, y2; 

{ 

double d2 = 9; 

double rl, r2: 

double sign; 

rl = (xl -x2)/2.0: 
r2 = (yl . y 2) / 2.0; 
if (rl > 0) sign = 1; 
else sign = -1; 
r2 = rPrl + r2*r2; 
r2 -=d2; 

if(r2<0)r2 = 0; 
•2 « sign * sqrt(r2); 

} 

/* **** minimal vector math ••••• •/ 

get_point(prompt, x, y) 
char 'prompt; 
double *x, *y; 
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{ 

printf(prompt); 
scanf("%If %lf\ x, y); 

} 

get_vector(prompt, v) 
char 'prompt; 
double *v; 

{ 

printf(prompt); 

scanf{"%lf %lf %lf\ &v[0], &v[l], &v[2]); 

} 

print_vector(prompt, v) 
char 'prompt; 
double *v; 

{ 

printf(prompt): 

printf("[%5.31f, %5.31f, %5.31f]\n", v[0], v[I], v[2]); 

} 

print_matrix(prompt, m) 
char 'prompt; 
double m[][3]; 

{ 

int i; 

for (i = 0; i < 3; r+-r) { 

prim_vector(prompt, m[i]); 

} 

} 

vadd(c, a, b) /• c = a + b V 

double *c, *a, *b; 

{ 

int i; 

for (i = 0; i < 3; i— ) c£i] = a[i] + b[i]; 

} 

vsub(c, a, b) /* c = a - b •/ 

double *c, *a, *b; 

{ 

int i; 

for (i = 0; i < 3; i++) c[i] - a[i] - b[i]; 

> 

vcmul(c, a, b) /• c * a * b */ 

App. A, Page 4 




WO 97/0674^H PCT/US96/13956 



double a; 
double •b, *c; 



} 



int i; 

for (i = 0; i < 3; H-+) c [i] = a • b[ij; 



vimul(c, a. b) I* c = a • b •/ 

double •a, »b. *c; 

10 { 

int i; 
•c = 0; 

for (i = 0; i < 3; i++) { 
•c += a[i] • b[i]; 

15 } 

J 



vemul(c, a. b) /* c = a X b */ 

double *c. 'a, »b; 



20 { 



double tmpl[3]: 
double tmp2[3]; 
int i; 



25 for (i = 0; i < 3; i++) { 

unplfi] = afi]; 
tmp2[i] = b[i]; 

) 

c[0] = tmpl[l] • tm P 2[2] - tmpl[2] • tmp2[l] 
30 c[l] = tmpl[2] • tmp2[0] - tmpl[0] • tmp2[2] 

c[2] = tmpl[0] • tmp2[l] - unpl[l] • tmp2[0] 

} 

vnorm(a) 
35 double a[3]; 



double r; 
int i; 



40 r - 0; 

for (i = 0; i < 3; i++) { 
r += a{i] • a[i]; 

) 

r = sqrt(r); 
45 for (i = 0; i < 3; { 

a[i] /= r, 

> 
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} 



vmmuKc. a. b) /* c = Ab ■/ 
^ double -c, a(][3J. «b; 

int i; 

double tmp(3]; 
for (i - 0; i < 3; i++) tmp(i] = bp]; 
vimul(&c[0], a[0], tmp); 
vimuK&cfl], a[l], tmp); 
vimul(&c[2], a[2], tmp); 



mmul(c, a, b) /• C = AB V 

^ double C Q[3], iOPJ. bQ[3]; 

int i, j, k; 

double tmpl(3J(3]: 

double tmp2(3][3J; 

for (i = 0: i < 3; i-^) { 
for (j = 0; j < 3: j++) { 

tmpimo] - ■[001; 
«mp2[i]D] = b[i][j]; 

) 
} 

for (i = 0: i < 3; i++) { 

for (j = 0; j < 3; j-*-) { 
c[i][j] = 0; 

for (k = 0; k < 3; k+-t-) { 

c[i]D] += «mpl[i][k] • tmp2[k][j]; 



} 



mtrans(a) 

double a Q(3]; 

{ 



double 



tmp; 



tmp = a[0][l]; a[0][l] = a[l][0]; a[l][0] = tmp; 
tmp = a[0][2]; a[0][2] = a[2][0]; a[2][0] = tmp: 
tmp = a[l][2]; a[l][2] = a[2][l]; a[2][l] = tmp; 
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mnorm(a) 

double aQ[3]; 

{ 

int i; 



for (i = 0; i < 3; { 
vnorm(a[i]): 

} 



dsl/206215.2 
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Appendix B 

c Brigham & Womens Hospital 1995 



rarc2.e 
T2 RARE 

K. Oshio 5-26-1993 



-5.2- 



5-26-1993 sc -> re 

7-5-1993 phase fix 

7-8-1993 rect. nopwTap 

15 7-8-1993 chem sat amp (flip_rfcssat = 110) 

7-9-1993 512. ysize fix 

7-10-1993 re-> rare2 

7- 27-1993 YPOS 
10-05-1993 <math.h> (SpSai etc) 

20 5.3 

1-18-1994 5.2 -> 5.3 (- indent) 

1- 18-1994 pass_lp(). slice _ip() 

2- 15-1994 48 eti (time_ssi - 4000) 
2-15-1994 a_gyl fix (eg_phaseres->yres) 

25 2-19-1994 watchdogcount increased (works) 

9-13-1994 3/4 fov etc for cholangiogaphy 

9-13-1994 rl-r9, rotate for plane correction 

2-22-1995 frame rotation 

30 4-4-1995 Prescan update (Prescan.e / grad_rf_data.h) 

4-4-1995 gzk fix 

4-8-1995 Prescan'gzk (works ... not sure about scaline) 

4-11-1995 scaling fixed (NtAX_ENTRY_POINTS « 9) 

8- 2-1995 Dixon 

35 8-5-1995 KOChemSAT 

plans 

bugs 



40 - nopwrap (7-15-95) 

- multi-pass ... ssi_time 

•/ 

45 @inline epic.h . 
©global 
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#include <maih.h> 
^include "em _psd_ermes. i n " 
^include "stddefh" 
^include M epicconf.h" 
5 ^include w pulsegen.h n 
#include Tiker.h" 

FILTER JNFO echoljilt: /• filter info stmct •/ 

10 /* PSD_FILTER_GEN psd_filt_spec[4] ; /• defined in epic.h •/ 

static char supfailfmtQ = "Support routine %s failed"; 

/* Array to hold max Bl for each entry point. V 
15 *Jefine MAX ENTRY_POINTS 9 

float maxB 1 [MAX_ENTRY_POINTS] ; 

float maxBlSeq; 

(ainline SpSat_blscaIe.e SpSatGIobal 
20 ^inline KOChemSat.e ChemSatGIobal 

int debugstate = 0; /• printdbgO on/off */ 

@ipgexport 

int offjrfcszfD AT A_ACQ_MAX] ; /* used in preprocessor •/ 

25 /• generated code •/ 

/• -CV- */ 
@cv 

(ginline SpSat_blscaIe.e SpSatCV 
30 @inline KOChemSat.e ChemSatCV 

/* main */ 



int time_ssi = 4000us: /• was 250us (orig), 2000 (5.2) •/ 

int queue_size = 4096; /• was 200 (orig). 1024 (5.2) •/ 

35 int watchdog_count = 3; /* (x5sec) default=l ? */ 

int tdO = G RAD_UP DA TE_T1 M E : /• cardiac delay •/ 

int tiead = 1ms; /• init dead time */ 

int pos_stan; /• start of pseq (excluding sat) •/ 

int dbg~= 1; /* same as 4x, 0.CPMG, 1 :RARE */ 

40 int use_ermcs: /• error message (should be in epich) V 

int bw_rfl; 

int bw_rf2; /• RF pulse bandwidth •/ 

int pw_ramp = 500usJ* default ramp width */ 

int pw_framp = 5 00 us; /* fast ramp */ 

45 int scanslot « 1; /• 8msec filter V 

float area_gzl ; 

float ob_target; 
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int trig_code = TRJG_LINE; 

im act_tr; 

int dda = 2; /• disdaq •/ 

int yres; /• actual YRES V 

5 int recno = -1; /• receiver -1 for all 0:L. 1:R */ 

int etc; /* effective TE •/ 

int nshot: 

int slicecm « 0; /• slices B4 pause (grass) •/ 

int slice_size: 

10 int dixon_de!ay = 0; 

/• CVs for Assymetric FOV •/ 

int transformsize; 

float phasefov; 



15 



20 



/• CVs for oddnex cases... •/ 

int oddnex: /• flag for oddnex case */ 



int 
int 



oddnex_npw; /• flag for oddnex and npw case */ 
evennexjipw; /* flag for evennex and npw case V 



/* rotation */ 

int rot_fix = 0; /* 0: no correction. I : correction */ 

float f_shift - 0; /• trans vector (mm) V 

float p_shift = 0; /* trans vector (mm) */ 

25 float z_shift = 0; /• trans vector (mm) •/ 

float rli - 1 .0: /• rot matrix •/ 

float r 22 = 1.0; /• rot matrix */ 

float r 33 = 1.0; /* rot matrix */ 

float r!2 « 0; /• rot matrix */ 

30 float r 23 = 0: /• rot matrix V 

float rl3 - 0: /• rot matrix •/ 

float r21 =0; /* rot matrix •/ 

float r 32 « 0; /* rot matrix •/ 

float r31 =0; /* rot matrix V 

35 



P rh cvs V 

shon autolock = 0; 

short blank = 4; 

short nograd = 0; 

40 short nofermi = 0; 

shon rawdata = 0; 

short saveinter = 0; 

short zchop = 1 ; 

shon eepf = 0; 

45 shon oepf = 0; 

shon eeff = 0; 

shon oeff « 0; 
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15 



20 



short nreceiver - 1; 

ini num_conc _grad = 3; 

/* adv panel •/ 
int min_seql; 
ini min_seq2; 
int min_seq3; 
int tmin_total; 
int max_gyl; 

/* prcscan •/ 

©inline Prescan.e PScvs 

short pre__pass « 0; /* ps pass* •/ 

shon pre_slice =*0; /* ps sliced •/ 

©host 

#include w sar_pm.h M 

^include "grad_rf_data.h w 

©inline Prescan.e PShost 

abstract(T2 weighted RARE"); 
psdname("rare2 H ); 



/• -CVINIT- V 
25 int 

cvinitQ 
{ 

#ifdef ERMES_DEBUG 
uscermes = 0; 

30 *else 

useermes = 1; 

#endif 

^include M cvinit.in" 

35 /• config var ? V 
EpicConf(); 

/• ink opcv •/ 

pititle « 0; 
40 piuset = useO; 

cvdesc(pitide, "User Variables"); 

opuserO = 0.0; 

cvmax(optr, 30s); 

cvmax(ihtr, 30s); 
45 cvmax(opshhick, 100); /• for projection V 

cvmax(opfast, 1 ); /* enable fast option •/ 

cvmax(opmph, 1); /• enable multi-phase option V 
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cvdef(opirmode. 0); /* non-sequential */ 

cvmax(oprbw. 32); P default max = 16 •/ 

oprbw = 32.0; 

cvdetfopetl, 7); 

cvdef(opte, 18ms); 

ete = opte • (opetl + 1) / 2; 

cvdef(opxres. 256); 

cvdef(opyres. 256); 

fh=l; /• init V 



/* init for oddnex flags •/ 
oddnex_npw = 0; 
oddnex = 0; 

evennex_npw » 0: 
15 tmenex = 0; 

/• ramp (-d = -a), others are set auto */ 

amptarget(&ob_targeL num_conc_grad): 
pw_gzrfla = pw_ramp: 
20 pw_gzrfl d = pw_ramp: 

pw_gzrGa = pw_ramc; 
pw_gzrf2d = pw_ramp; 

/* trans spoiler */ 
25 pw_gzk = 4ms; 

pw_gzka = pw_ramp: 
pw_gzkd = pw_ramp; 

/• RF (has to be consistent with gradjf.h) */ 
30 /• (these have to be initialized even if not used */ 

pwjfl = 3200; 

pw_rf2 = 3200; 

pw_rf22 = 3200; 

pwjgzrf 1 = pw_rfl ; 
35 pw_gzrf2 = pw_rf2; 

pw_gzrf22 = pw_rf2: 

cycjfl = 1.0; 

cyc_rG = 1.0; 

cyc_rf22 = 1.0; 
40 gscale_rfl = 0.9; 

gscale_rf2 = 0.9; 

gscaie_rf22 = 0.9; 

/• readout •/ 

45 echol_fiIt = f!6; /• copy struct V 

/* advisory panel */ 
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piadvise = 1; 

piadvmin = (1 « PSD_ADVTE) + (I « PSD_ADVTR) + (1 « PSD_ADVFOV) 
piadvmax = (1 « PSD_ADVTE) + (1 « PSD_ADVTR) + (1 « PSD_ADVFOV) 
piadviime = (1 « PSD_ADVMINTSCAN) + (1 « 
PSD_ADVMAXLOCSPERACQ) 

+ (1 « PSD_ADVMINACQS) + (1 « PSD_AJDVMAXYRES); 
avminslquant = 1; 
avmaxte = 500ms: 
avmaxtr = 10s; 
avmaxfov = 480; 

/• SAT V 

vrgsat = SINC_SAT; 
SpSatlnit(vrgsat);* 

/• ChemSAT •/ 

@inline KOChemSat.e ChemSatlnit 

/* PS •/ 

PScvinitQ; 
AScvinit(); 

return SUCCESS: 
} /• cvinit() •/ 

©inline SpSat_blscale.e SpSatlnit 
©inline SpSat_blscale.e SpSatCheck 

/* -FILTER- •/ 
int 

calc_fiher(info, bw\ res. slot)/* make filter (move to lib) V 
FILTER JNFO *info; /• output V 
/• bandwidth •/ 
/• xres V 



float 


bw: 


int 


res; 


int 


slot; i 


float 


abw; 


int 


taps; 


int 


prefill; 



PSD_FILTER_GEN 'spec; 

/* calc info •/ 

if (bw <= 0.0) { 

info->fslot = 0; 
return FAILURE; 
} else { 

info->decimation = 64 / (int)bw; 
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} 



if (info->decimation = 0) 

info->decimation - I; 
abw = 64.0 / (float)info->decimation; 
info->bw = abw • 10.0 / 10.0; /* round to 0.1 kHz */ 
info->tdaq = res * info->decimation * PSDJTSP; 
info->fslot = slot: 
info->fname = 0; 



10 /• filter gen V 

if (info->fslot < 4 || info->fslot > 8) 
return FAILURE; 

spec = &psd_fik_spec[0]; 

psd_filt_spec[l].filter_slot =0; 
IS psd_fiIt_spec[2].fiiter_slot = 0; 

psd_filt_spec[3].filter_slot = 0; 

switch (info->decimation) { /• Empirical measured values */ 
case 1: 

taps = 67; /• not tested yet •/ 
20 break; 
case 2: 

taps = 67; 

break; 
case 3: 

25 taps = 102; 

break; 
case 4: 

taps = 180; 

break: 
30 case 5: 

taps = 201: 

break; 

default: /■ Anything past 255 taps. Who cares? */ 

taps = 255; 
35 break; 
} 

prefill « (taps + 1) 12 - 1; 

spec->filter_slot - info->fslot; 
40 spec->taps = taps; 

spec->dec = info->decimation; 

spec->outputs = res; 

spec->prefills = prefill; 

spec->flush = 0; 
45 spec->tflush = 0; 

spec->post_flush = 0; 

spec->queuing = 0; 



App. B. Page 7 



WO 97/06744 



PCT/US96/13956 



spcc->rshift = 0; 
spec->alp = 0.46; 

if (spec->dec & 0x01) /• quadrature detection •/ 

spec->det = 0; /• if odd must do hardware detection •/ 

else 

spec->det = 1; /* if even can do coefficient detection •/ 
if (spec->dec = 1) /* for filters with no decimation •/ 

spec->dfg = 1.0/* must use unitary gain */ 
else /• otherwise use gain of 2.0 •/ 

spec->dfg = 2.0;/* This is a bug that will be fixed. */ 
f* (says Bruce Collide) V 
spec->filfrq = bw; 
spec->picw = 0.985; 
spec->wt = 1.0: 
spec->invert_q = 0; 
return SUCCESS; 
} /* calc_filter() •/ 

/* -CVEVAL- •/ 
int 

cveval() 
{ 

int status: 

int av_temp_int; 

/• screen control */ 
/• flip angle •/ 

pifanub = 0; 

pifavatt = 10; 

pifava!3 = 20; 

pifaval4 = 30; 

pifavalS = 45; 

pifava!6 = 60; 

/• NECHO •/ 

piechnub = 0; 
piechvall = 1; 
piechval2 = 2; 

/• TE •/ 

pile 1 nub = 63; 
pitelval2 = 12ms; 
pitelva!3 = 15ms; 
pitelval4 = 18ms; 
pitelvalS = 20ms; 
pitelval6 = 30ms; 
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/• ETL •/ 

pietlnub 5 

pieUval2 

pieilva!3 

pietlval4 

pietlvalS 

pietival6 

/• TE2 V 

pite2nub 



63; 
= 5: 
: 7; 
s 11; 
■ 15: 
= 64; 



= 0; 



/* TR */ 

piimub = 6; * /• number */ 
pitrvaI2 = 500ms; 
pitr\aI3 = 2000ms: 
pilrval4 = 2500ms 
pitrvalS = 3000ms 
pitrva!6 = 4000ms: 

r NEX V 

pinexnub =57; /•#..###•/ 
pinexva!4 = 1; 
pinexva!5 = 2; 
pinexval6 = 4; 

/* XRES/YRES V 

pixresnub = 6; /*.£....*/ 
piyrcsnub =62; 

/• bandwidth */ 

pircbnub = 0; /* number */ 

pircb2nub = 0; 
pircbval2 = 32.0; 
pircbva!3 = 16.0; 

/* multi-phase •/ 

if (opmph = PSD_ON) { . 

pimphscrn = 1; 

pifphasenub = 6 j* number */ 

pifphaseval2 = 1; 

pif^)haseva!3 = 2; 

pifphaseval4 = 5; 

pifphasevaI5 = 10; 

pifphaseval6 = 15; 
} else { 

pimphscrn = 0; 

pifphasenub = 0; 
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} 

if (opfphases = I ) { 

pisldeinub = 0; 

piacqnub = 0; 
} else { 

pisldeinub = 6; 
pisldelval3 ■ 500ms; 
pisldelval4 = 1000ms; 
pisldelval5 = 2000ms; 
pisldelval6 = 5000ms; 

piacqnub = 2; 

) 

/• Variable FOV buttons on or off depending on square pixels */ 

if ((exist(opsquare) = PSDJDN) || (exist(opnopwrap) = PSD_ON)) { 

piphasfovnub = 0; 
} else { 

piphasfovnub = 7; 

} 

/**•• Asymmetric Fov ••••/ 

/* handling for phase (y) resolution and recon scale factor.*/ 

if ((exist(opnopwrap) = PSD_ON) existcv(opnopwrap)) { 

rhphasescale = !.0: 

setexist(opphasefov, PSD_ON); 

_opphasefov.fixedflag = 0; 

opphasefov = 1.0: 

_opphasefov.fixedflag = 1; 

eg_phaseres = exist(opyres); 
} else if ((exist(opsquare) = PSDJDN) existcv(opsquare)) 

{ 

rhphasescale = (float)exist(opyres)/( float )exist(opxres); 

setexist(opphasefov,PSD_ON); 

_opphasefo\\fixedflag = 0: 

opphasefov = rhphasescale; 

_opphasefov.fixedflag = 1; 

eg_phaseres = exist(opxres); 

} 

else if ((exist(opphasefov) != 1.0) 

&& existcv(opphasefov) 
&& (exist(opsquare) != PSDJDN) 
&& existcv(opsquare)) { 
rhphasescale = exist(opphasefov); 
eg_phaseres = exist(opyres); 
} else { 
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rhphasescale = 1.0; 
_opphasefov.fixedflag = 0; 
opphasefov = 1.0; 
_opphasefov.fixedflag = I; 
eg_phaseres = exist(opvres); 

} 

/• multi-coil •/ 

if (recno = -1) { 

rhdabOs = cfrecvst: 

rhdabOe = cfrecvend; 
} else { 

rhdabOs = recno; 

rhdabOe = recno; 

} 

nreceiver = rhdabOe - rhdabOs + 1; 

/* slice order type */ 

if (seqtypef&seqjype) = FAILURE) { 
epic_serror( "seq t ype( ) fai 1 ed " ) ; 
return FAILURE; 

} 

acq_rype - TYPSPIN; 

/• chem sat •/ 

if (ChemSatEval(&cs_sartirne) — FAILURE) { 
epic_serrorC rt ChemSatEval() failed"); 
return FAILURE; 

} 

/* KOChemSat.e will handle these 

cs_sanime = RUP_GRD(pwjfcssat + pw_gykcsa + pw_gykcs + pwjgykcsd) 

+ 1ms; 
if (!cs_sat) 

cs_sanime = 0; 
flip_rfcssat = 110.0; 



/• sp sat •/ 

if (SpSatEvaI(&sp_sattime) = FAILURE) { 
epic_serror( M Spatial Sat eval failed"); 
return FAILURE; 

} 

/• RFU RF2 V 

bw_rfl = 4.0 • cyc_rfl / (pwjfl / (float) I. Os); 
bw_rf2 = 4.0 • cyc_rf2 / (pw_r£ / (float) 1.0s); 
if (ampsHce(&a_gzrfl % bw_rfl, opslthick. 

gscale_rfl, TYPDEF) = FAILURE) { 
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epic_serror( M ampslice failure for gzrfT); 

if (ampslice(&a_gzrf2. bw_rf2, opslihick. 

gscale_rf2, TYPDEF) = FAILURE) { 
5 epic^scrrorC'ampslice failure for g2rf2 n ); 

a_gzr£22 « a_gzrf2: 

/* crusher / ramp V 
10 /♦ optramp(&pw_gzrfla, a_ezrfl, 1, TYPDEF); V 
pw_gzrfld = pw_gzrfla: 

pw_gzrf2r = 1200; 
pw_gzr£221 = 1200; 
15 pw_gzr£22r = 1200: 

a __gzrf2r = 0.9; 
a_gzr£221 = 0.9; 
a_gzrf22r = 0.9; 

optramp(&p\v_gzrf2ra. a_gzr£2r. 1, TYPDEF); 
pw_gzrf2rd = pw_gzrf2ra; 



20 



25 



30 



pw_gzrf2a = pw_gzrf2ra; 
pw_gzrf2d = pw_gzrf2ra: 
pw_gzrf21a = pw_gzr£2ra; 
pwjgzrf21d = pw_gzrf2ra: 
pw_£zrf22a = pw_gzrf2ra: 
pw_gzrf22d = pw_gzrf2ra: 
pwjgzrG21a = pw_gzr£2ra: 
pw_gzrf221d = pw_gzrf2ra; 

pw_gzrf22ra = pw_gzr!2ra: 
pw __gzrf22rd = pw_gzrf2ra; 



area_gzl = (pw_gzrfl + pw_gzrfld) • a_gzrfl / 2.0; 
35 if (amppwlcrsh(&gradz[GZRF2L_SLOTj, &gradz[GZRF2R SLOT], 

area_gzl, a_gzrf2. objarget) = FAILURE) { ~ 
epic_serrorCainppwlcrsh failure for gzrG"); 

/• read bandwidth / filter •/ 
40 if (exist(opfov) < 160) 

oprbw = 16.0; 

else 

oprbw = 32.0; 

if (oprbw = 16.0 &.&. exist(opxres) = 256) { /• default filter •/ 
45 scanslot = 1; 

echoljilt = fl6; 
} else { /• real time fill gen •/ 
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scanslot = 4; 

calc_fiIter(&echol_filL oprbw. exist(opxres), scanslot); 

) 

5 /• read/phase amp */ 

if (optr > 3s) dda = 0; 
dsc dda = 2: 

if (opetl > 32) timejsi = 8ms: 
10 else time_ssi = 4ms; 

/* nop •/ 

if (exist(opnopwrap) = PSDJDN) nop = 2; 
else nop = 1 ; 

15 

/• yres / nshot •/ 

rhnframes = opyres * fn * nop * opphasefov; 
rhnframes = rhnframes / opetl / 2 * opetl * 2; 
if (rhnframes = 0) rhnframes - opetl; /* single shot •/ 
20 yres = rhnframes / opphasefov; 

if (dbg — 0) { /• CPMG •/ 

nshot = rhnframes: 

rhnecho = opetl; 
25 ete = exist(opte); 

} else { /• RARE V 

nshot = rhnframes / opetl; 

rhnecho = opnecho: 

ete = exist(opte) * (opetl + 1) / 2; 

30 } 

if (amppwencode(&a_gy 1 . &pw_gyl. /* nop is float •/ 
exist(opfov), yres / (int)nop. 16ms, 3) = FAILURE) { 
epic_serror("amppwencode call failed"); 

35 ) 

if fampfov(&a_gx\v. oprbw. exist(opfov)) = FAILURE) { 
epic_serror("ampfov call failed"); 

} 

pwjgxw = echol_filt.tdaq; 
40 a_gyla = ajgyl; 

a_gyl = -a_gyla; 
pw_gyla = pwjgyl; 

if (optramp(&pw_gxwa. a_gxw, 1, TYPDEF) = FAILURE) { 
epic_serror( M optramp call failed"); 

45 , 

pw_gxwd = pw_gxwa; 



App. B. Page 13 



WO 97/067< 




PCT/US96/I3956 



/* pause (grass) •/ 

if (exisi(opslicecni) == 0) { 

slicecnt = exist(opslquani); 
} else { 

slicecnt = exist(opslicecnt): 

} 

/* adv panel */ 

avminnecho = 1; 
avmaxnecho = 1; 

min_seql = pwjfl / 2 + pw_ezrfld /* 90 - 180 •/ 
+ pw_gzrf2la + pw_gzrf21 
+ pw_gzrf2a + pw_gzrf2 / 2; 
min_seq2 = pw_rf2 + pw_gyl * 2 + pw_g X w; /• 180 - 180 */ 
min_seq3 = pw_gzrf2 / 2 .'* 180 - read 

+ pw_gzrf2ra + pw_gzrGr + pw_gzrf2rd 
+ pw_gxw / 2; 
if (mintel(&avminte, eg_phaseres, exist(opfov) % 
min_seql, min_seq2. min_seq3 % 
TYPSPIR TYPNVEMP. echol_fiiudaq. 
pw_rfl, pw_rf2. TYPNFC 1.0) = FAILURE) { 
epic_serror("mintel failed"); 
return FAILURE; 

} 

avmaxte = eie; 

tmin_total = tlead + tdO sp_sattime - c$_sattime 

+ pw_gzrf 1 a + pw_rfl / 2 

+ exist(opte) * opetl 

+ pw_gxw / 2 + pw_gxwd 

+ pw_gzka + pw_gzk pw_gzkd: 
if (mintr(&avmintr. seq_type. 

(long)tmin_totaI. 1. trig_code) = FAILURE) { 

return FAILURE: 

} 

max _gyl = exist(opte) - pw_rfl / 2 - p\v_gxw / 2; 

if (minfov(&avminfov, exist(opfov), seq_jype, max_gvl % Is, 

echol Jilt.bw, exist(opyres), TYPNFC p\v gxw. L0) = FAILURE) 
return FAILURE; 

} 

if (maxyres(&av_tempjnu ob_target, max_gvl) = FAILURE) { 
return FAILURE; 

} 

if ((avmaxslquam = (int)(optr / avmintr)) = 0) 

return FAILURE; 
slice_size = (1 + rhnframes rhhnover) * 2 

* rhptsize • rhfrsize * exist(opnecho) * nreceiver; 
maxslquanttps(&av_iemp_int, 
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(int)rhimsize. 1, (int)rhdayres. slice_size); 
if (av_temp_int < avmaxslquant) 

avmaxslquant = av_temp_int: 
if (opirmode = PSD_SEQMODE_ON) 

avmaxslquant = 1; 
acqs = (opslquant - 1) / avmaxslquant + 1; /• # of passes 
sliceinl(&slquantl, acqs. acq_type); /• # of slices in pass 1 
avmaxacqs = acqs; 

avmintscan « optr * nshot; /* time for 1NEX V 
pisctim4 = avmintscan; 
pisctimS = avmintscan • 2; 
pisctim6 = avmintscan * 4; 

/* check current value */ 

if (exist(optr) < avmintr jj exist(opte) < avminte) 
return FAILURE; 

/* prescan •/ 

if (PScveval() = FAILURE) return FAILURE: 
if (AScvevalQ = FAILURE) return FAILURE; 

return SUCCESS; 
} /• cveva]() V 

/• -CVCHECK- •/ 
int 

cvcheck() 
{ 

if (exist(opte) < avminte 

|| exist(optr) < avmintr 

|| exist(opfov) < avminfov) 
/• this shouldn't happen, checked in cvevalQ */ 

return FAILURE; 

return SUCCESS; 
} /* cvcheck() V 

r -PREDOWNLOAD- V 
int 

prcdownload() 
{ 

int i; 

sp_satstan = RUPJ3RD((int)tIead + GRAD_UPDATE_TIME); 
cssatstan = sp_satstart + sp_sattime + 1ms: 
pos_stan = sp_satstart + sp_sanime + cs_sartime: 
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/• multi-slice •/ 

if (orderslice(seq_rype. opslquanu 

slquantl. trig_code) = FAILURE) { 
epic_serror( H orderslice call failed"); 

} 

pitr = psl_tr; /♦ PS1 TR •/ 

pichop = 0; /* No chop •/ 

pitspl = Q.decimation * PSD_TSP; /• Sample lime for CFL •/ 
pitsp2 = f500.decimation • PSDJSP; /• Sample time for CFH V 
pislquant = slquantl; /• U of 2nd pass slices •/ 

baseline = 0; 
nex « opnex; • 
exnex = opnex: 

/* scan time •/ 

act_tr = RDN_GRD(optr / slquantl); 
if (opslicecnt = 0) { 

pidmode = PSD_CLOCK_NORM: 
} else { 

pidmode = PSD_CLOCK_PAUSE; 

) 

pitscan = (nshot + dda) * optr * opnex * acqs; /• clock */ 

/* image anotation V 
ihtr = optr; 
ihtel = ete; 
ihflip * opflip: 
ihnex = opnex: 
ihvbwl = oprbw; 

/* entry table •/ 

if (entrytabinit(entrv_poim table. (int)ENTRY_POINT_MAX) 
= FAILURE) { 

epic_serror("Can*t initialize entry point table."); 
return FAILURE; 

} 

entry joint_table[L_SCAN].epfiIter = scanslot; 
entry_poim_table[L_APSl].epfiher = F16_SLOT; 
entry joint jable[L_CFH].epfilter = F500_SLOT; 
entryjoint_tabIe[L_CFL].epfilter = F2_SLOT; 

entry_point_table[L_MPS 1 ] = entry_point_table[L_APSl]; 
entry_point_tabIe[L_APS2] = entry _poim_table[L_MPS2] = 

emry_point_table[L_SCAN]; /• copy scan into APS2 & MPS 2 •/ 
strcpy(entry_point_table[L_SCAN].epname. "scan"); 
strcpy(entry_point_table[L_CFH].epname, "cfh"); 
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stirpy(enir>'_ooini_tab!e[L_CFL].epname. "cfT): 
strcpy(eniry'_poini_iabIe[L_APSl].epname. "apsF): 
strcpy(entr\'_poim_iabie[L^MPSl ].epname. "mpsP); 
strcpy(emryjoim_tabie[L_APS2].epname. *aps2"); 
strcpyCenuyjDoim^tablefL^MPSZJ.epname. ,, mps2 ,,, ); 

^include "predownioad.in" 

/• !!! NOT included in macro !!!! •/ 
ia_gylw = ia_jy 1 ; 
ia_gylaw = ia_gyla: 

/• RF pulses are scaled by entry point basis */ 
maxBlSeq = 0.0: 

for (i = 0; i < MAX_ENTRYJ>OINTS; i++) { 

if (peakfll(&maxBl[i] ? i. RF_FREE. rfpulse) = FAILURE) { 
epic_error(use_ermes. supfailfmt. 

EM_PSD_SUPPORT_FAILURE 
1. STRING_ARG, H peakBP) 
return FAILURE; 

} 

if (maxBl[i] > maxBlSeq) 

maxBlSeq = maxBl[i]; 

if (setScale(L_SCAN\ RF_FREE. rfpulse. maxBl[L SCAN]. 

maxBl[L_SCAN] / maxBlSeq) = FAILURE) { 
epic_error(use_ermes. supfailfmt, EM_PSD_SUPPORT FAILURE 

1, STRING_ARG. "setScale'T 
return FAILURE; 

} 

ia_rfl = maxjDgJamp • Crfpulse[RFI_SLOT].amp): 
ia_rf2 = max_pg_iamp • CrfpulsefRF2_SLOTJ.amp): 

/• debug 

if (exist(opsatx) > 0) 

iajfsxl = max_pg_iamp • CrfpulsefRFSXl SLOT].amp); 

ia_rfsx2 = maxjjgjamp • Crr P ulse[RFSX2^SLOT].amp)- 
if (exist(opsaty) > 0) 

ia_rfsy! = max_pg_iamp • CrfpuisefRFSYI SLOT].amp); 

ia_rfsy2 = maxjgjamp • CrfpuIse[RFSY2^SLOT].amp);' 
if (exisi(opsatz) > 0) 

iajfszl = max_pgjamp • (*rfpulse[RFSZl_SLOT].amp); 

iajfsz2 = max_pg_iamp • (Mpulse[RFSZ2_SLOT].amp);' 



/• sat •/ 

SpSaUAmp(O); 
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SatPlaccmcnt(acqs); 
if (cs_sat == PSD_ON) 

iajfcssat = max_pg_iamp • (Mpulse[RFCSSAT_SLOT].amp); 

/* prescan •/ 

if (prescanslice(&pre_pass. &pre_slice, opsiquant) = FAILURE) { 
cpic - serror( ,, Can't select a slice for prescan."); 
return FAILURE: 

} 

PSpredownload(); 
ASpredownloadQ; 

/* ===== recon header variables ============ */ 

/* rhrype •/ 

rhtype = 0; 
if (nex > 1) 

rhtype j= RHTYPCHP; 
if (opimode — PSD_CINE) 

rhtype |= RHTYPCINE; 
if (oppseq — PSDGE || 
oppseq = PSDJTOF || 
oppseq = PSD_TOFSP || 
oppseq = PSD_PC || 
oppseq — PSDJPCSP) 
rhtype |= RHTYPGR; 
if(fn = 0.5) 

rhtype |= RHTYPFRACTNEX; 
if (oprect) 

rhtype |= RHTYPSTRIP; 
if (opimode = PSD_3D) 

rhtype |= RHTYP3D; 
if (nop = 2 && oppomp = PSD_OFF) 

rhtype |= RHTYPNPW; 
if (pitfeextra > 0) 

rhtype |= RHTYPFRACTECHO; 
if (fh — 0.75) 

rhtype |= RHTYP75NEX; 
if (oppomp) 

rhtype |= RHTYPPOMP; 

/• rhformat •/ 

rhformat = nograd + 2 • nofermi; 
if (opimode = PSD_3D) { 
if (zchop) 

rhformat |= 4; 

else 

rhformat |= 8; 
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} 

if (opscic — PSD_ON) 
rhformat |= 16; 



/• others •/ 

rhbline = baseline; 
rhblank = blank; 

rhfrsize = exist(opxres); /* rhnframes is in eval •/ 

rhnslices = opslquant * opfphases; 

rhptsize = opptsize; 

rhnavs = lMax(2, nex / 2, 1); 

rhdacqctrl = rawdata + 

2 * eepf + /• even echo phase flip •/ 
4 • oepf + /* odd echo phase flip */ 
8 * eeff + /• even echo freq flip •/ 
16 * oeff; /* odd echo freq flip •/ 

rhexecctrl « 9 + 2 * autolock + 64 • saveimer; 

rhvquani = opvquant; /• VQUANT ... does this work ? •/ 

rhslblank = pislblank; 

rhzeroph = (yres • nop / 2) / (I + opreci) + 0.5; 
if (pitfeexira > 0) 
rhnwin = 8 

else 

rhnwin = 0 

if (pitfeexira > 0 || (fn — 0.5) && (nop = 2)) 
rhntran = 4 

else 

rhntran = 2 



rhfermr - rhfrsize / 2; 
rhnpasses = acqs; 



rhrawsize « slquantl * slice_size; 

rhrcyres = exist(op\Tes); 
rhrcxres = exist(opxres); 
rhmethod - oppomp; 
rhdaxres = exist(opxres); 
rhdayres = rhnframes + rhhnover + 1; 
if (exist(opxres) > 256 || exist(opyres) > 256) 
rhimsize =512; 

else 

rhimsize = 256; 
rhnpasses = acqs; 

/* calculation for rhyoff with asym fov •/ 
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phasefov = exist(opphasefov) • exisr(opxres); 
transformsize = I; 

while (transformsize < phasefov - 0.001) transformsize «= I; 

if (transformsize > phasefov + 0.001) transformsize = exist(opxres); 

return SUCCESS: 
} /* predownload() ■/ 

©inline KOChemSat.e ChemSatEval 

©inline KOChemSat.e ChemSatCheck 
©inline SpSat_blscale.e SpSatEval 
©inline SpSat_blscale.e SatPlacement 

/* -PULSEGEN- •/ 
©rsp 

/* used by pulsegen & rsp •/ 

CHAR • entry_riameJist[ENTRY POINT MAX] = { 

"scan", 
"cfl", 
M cfh'\ 
"mpsl\ 
"apsF, 
n mps2". 
"aps2'\ 
"autoshim". 
0}; 

mt *rfl_freq. *rf2_freq; 

float •rfljhase. *rf2_phase; 

int *receive_freql; 

int cs_satindex « 0;/* sat index */ 

int sp_satindex = 0;/* sat index */ 

int deadtime; /• dead time for scanQ •/ 

int seqtime: /• time/seq */ 

short rsprot_orig[DATA_ACQ_MAX][9]; /• rotation matrix for SAT •/ 

/• array to store calculated phase offset and sign •/ 
typedef struct { 

int ysign; /• sign of the phase offset •/ 

int yoffs; /• calculated phase offset from rsp info structure */ 
} PHASEJ5FF; 

PHASEJ3FF phase_oflIDATA.ACQ.MAX]; 

int yres_phase; /* offset in phase direction in mm */ 

int yoffs 1; /• intermediate phase offset variable •/ 

WF_PULSE •echol; 

float orig_rot[3][3]; /• original rotation matrix */ 
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RSPJNFO •rsp_shift; /• copy rsp_info. and shift •/ 
int 

ssisatQ 
{ 

int next_slice; 

next_slice = sp_sat_index; 

sp^update^rot^matrixC&rsprot^origfnext^siicejlO], sat_rot_matrices, 
sat_rot_ex_num. sat_rot_df_num); 

/• if (opexor = PSD_ON) exorcist_ssi(); •/ 

return SUCCESS; 

} 

int 

dummyssi() 
{ 

return SUCCESS; 

J 

@Pg 

puisegen() 
{ 

int echo: /* echo index */ 

int rf2Ioc; /• center of 180 •/ 

char pname[256]; 
int i, j t ix; 

/* set watchdog timer */ 

serwatchdogrsp( watchdog_count); 

/* trigger delay (modify period to put delay) V 

WAIT(XGRAD, x_td0. Head. GRAD_UPDATE_TIME); 
WAIT(YGRAD. y_td0, tlead. GRADUPDATE TIME); 
WAIT(ZGRAD, z_td0. tlead. GRAD_UPDATE~TIME); 
WAIT(RHO, rhojdO. tlead. GRAD UPDATE JT1ME); ' 
WAITfTHETA, thetajdO, tlead, GRAD_UPD ATEJTI ME) ; 
WAIT(SSP, ssp_td0, dead GRAD_UPDATE_T1ME); 

r SpSAT •/ 

sp_satindex = 0; 

SpSatPG(vrgsat, sp_satstan. &sp_satindex, sp_satcard_loc); 

/• ChemSAT V 

cs_satindex = 0; 
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if (cs_sat) ChemSatPG(cs_satstart. &cs_satindex): 



/• 90 / dephaser V 
SLICESEL2(rfl, 

pos_stan + pw_gzrfla. 

pw_rfl. opslthick. opflip, cyc_rfl, TYPNDEF); 
TRAPEZOID(XGRAD, gxl, 



a_gxw • (pw_gxw + pw_gxwa) / 2.0. TYPDEF); 



/* alloc mem for DAB packets */ 

echol = (WF_PULSE *) A!JocNode((opetl + 2) * sizeof(WF_PULSE)); 

/• echo train •/ 

rf21oc = pmid(&rfl, "gzrfl", 0) + opte / 2; 
for (echo = 0; echo < opetl; echo++) { 
if (echo = 0) {/• first echo V 
SLICESELZ(rf2. 

rf21oc - pw_gzrf2 / 2. pwrO. 



opslthick. 180. cyc_rf2. TYPNDEF); 
TRAPEZOID(ZGRADB. gzrGl. 

pbee(&rf2. "gzrf2a". 0) - pw_gzrf21. 

0. TYPNDEF); 
TRAPEZOID(ZGRADB. gzrfir. 

pendf&rf2. "gzrf2d". 0), 

0, TYPNDEF); 



SLICESELZ(rf22. 

rf21oc - pw_gzrf2 / 2. pw_rf2. 

opslthick. 180. cyc_rf2. TYPNDEF): 
TRAPEZ01D(ZGRADB. gzrf221. 

pbee(<trf22. "gzrf22a". echo - 1) - pw_ezr021. 

6. TYPNDEF); 
TRAPEZOID(ZGRADB, gzrf22r. 

pend(&rf22. "gzrf22d-. echo - 1). 

0, TYPNDEF); 



TRAPEZOID(XGRAD. gxw. 

rf21oc + opte 12- pw_gxw 12 + dixon_delay, 

0, TYPNDEF); 
SlNUSOID(YGRAD. gyl, 

pbeg(&gxw. "gxw". echo) - pw _gyl, pw_gyl, 

ajgyl, 6, 0.0, 0.5. 0.0); 
SINUSOID(YGRAD, gyl a. 

pend(&gxw, "gxw", echo). pw_gyla. 



/* 



seqcore 



pend(&rfl, "gzrfl", 0) + pw_gxla. 



} else { 
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a_gyla. 0. 0.0, 0.5, 0.0); 



/• DATA ACQ V 

sprintf(pname, "echoMyMEd", echo); 

pulsename(&(echol[echo]), pname); 

acqq(&(echol(echo]), 

pbeg(&gxw, w gxw'\ echo) + psd_grad_acq_dela\\ 
scansloi. recno, DABNORM, PASSTHROUGH); 

rf21oc opte; 



/• killer V 

TRAPEZOID(ZGRAD, gzk, 

pend(&gxwd. "gxwd", opetl - 1) + pw_gzka. 
0, TYPNDEF); 

/• define sequence •/ 

SEQLENGTH(seqcore. act_tr ? seqcore): 
getperiod(&deadtime, &seqcore. 0); 
seqtime = act_tr - deadtime; 

/• prescan */ 

PSpulsegen(); 
ASpulsegen(); 

/* pass •/ 

PASSPACK(pass_pulse ? 49ms); /♦ tell Signa svstem we're done V 
SEQLENGTH(seqpass, 50ms, seqpass); 

/•sequencer memory alloc (original size) •/ 



ipg_alloc_instr[0] 


= 2048 


/• 


XGRAD (816) V 


ipg_alloc_insir[l] 


= 2048 


/• 


YGRAD (816) V 


ipg_alIoc_instr{2] 


= 2048 


/* 


ZGRAD (816) V 


ipg_alloc_instr[3] 


= 2048 


/• 


RHOl (816) V 


ipg_aIloc_instr[4] 


= 2048 


/* 


RH02 ( 64) •/ 


ipg_alloc_instr[5] 


= 2048 


/• 


THETA (816) •/ 


ipg_alloc_instr[6] 


« 2048; 


/• 


OMEGA (816) •/ 


ipg_alJoc_instr(7] 


= 2048; 


f 


SSP (1024) •/ 


ipg_alIoc_instr[8] 


= 64; 


/* AUX ( 64) V 



/• build instruction •/ 

if (buildinstr() = FAILURE) 
return FAILURE; 

/• Calculate the RF & slice frequencies •/ 

rfl_freq = (int *)AilocNode((opslquant * 2) • sizeof(int)); 

rf2_freq = (int *)AllocNode((opslquant + 2) * sizeof(int)); 
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rfl _phase = (float *)AllocNode( (opslquant + 2) • sizeof( float)); 
rf2_phase * (float •)AlIocNode((opsiquant + 2) * sizeof(float)); 
receivejreql = (int *)AIIocNode((opslquam + 2) * sizeof(im)); 

/* rotation / translation (modify rsp_info before setting up freq etc.) */ 
/* save original rot matrix (rsprot will be overwrinen) •/ 
for (i = 0; i < 3; i++) { 

for(j = 0;j<3;j++) { 
ix - j + 3 • i ; 

origjrot[i][j] - (float)rsprot[0][ix]; 

} 

} 

/• shift rspjnfo (x, y t z) •/ 

rsp_shift = (RSP INFO # )AllocNode(opslquant • sizcoffint)): 
if (rotjix) { 

for (i = 0; i < opslquant: { 

rsp_shift[i].rsptloc = rsp_info[i].rsptloc + z_shift; 
rsp_shift[i].rsprloc = rsp_info{i].rsprloc + f_shift: 
rsp_shift[i].rspphasoff = rsp_info[i].rspphasoff + p_shift; 

) 

} else { 

for (i = 0; i < opslquant; i++) { 

rsp_shift[i].rsptloc = rsp_info[i].rsptloc; 
rsp_shift[i].rsprloc = rsp_info[i].rsprioc; 
rsp shiftfij.rspphasoff = rsp_info[i].rspphasoff; 

} 

} 

/* Set the Slice Frequency */ 

setupslices(rfl_freq. rsp_shift. opslquant, 

a_gzrfL (float) 1. TYPTRANSMIT); 
setupslices(r£2_freq, rsp_shift. opslquant. 

a_gzrf2. (float) K TYPTRANSMIT); 
/* multi-slice CPMG phase V 

set_cpmg_phase(rf]_phase. rfl_freq, 0.0, pw_rfi / 2 + 9); 
set_cpmg_phase(rf2_phase. rf2_freq, PI / 2.0. pwjf2 / 2 + 9); 

seiupslices(receive_freql, rsp_shifi opslquant. 0.0, oprbw. TYPREC); 
setupphasetable(viewtable, TYPNORM, rhnframes); /• replace this •/ 

return SUCCESS; 
} /* pulsegen() V 

void 

set cpmg phascfphase. freq, offs. del) 

float *phase; /* phase array (output) */ 

int •freq; /• freq array (input) •/ 
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float offs; /* phase offset •/ 

int del; /* time from reset to center of pulse •/ 

{ 

int i; 

5 float ph; 

/* doesn't matter as long as same freq are used for rfl & rC */ 
/• necessary for innner volume, though */ 

1° for (i = 0; i < opslquant: i++) { 

ph » del • freqfi] / 0.59604648; /• phase •/ 
ph -= (int)ph; 

phasefi] = -ph • 2.0 • PI + offs; 

} 

15 } 

©inline KOChemSat.e ChemSatPG 
©inline SpSat_blscale.c SpSatPG 
©inline Prescan.e PSipe 
20 r End of ©pg V 

/• -RSPVAR- •/ 
@rspvar 

©inline Prescan.e PSrspvar 
25 ©inline KOChemSat.e ChemSatRspVar 
©inline SpSat_blscale.e SpSatRspVar 

shon viewiable[513]; /• max view t l (base 1) V 

shon aco_pir[256]; 
30 shon slc_in_acq(256]; 

int dabop; /• dab op (load, add etc) */ 

int acq_echol; /• dab on/off •/ 

int rspenu rspdda, rspbas. rspvus. rspgyl, rspasi: 

int rspesl, rspchp. rspnex. rspslq, rspsct; 

35 int rfl Phase, seqCount; /• SPGR V 

shon chopamp; 

/• -RSP- V 
©rsp 

40 ©inline KOChemSate CsSatMod 

©inline SpSat_blscaie.e SpSatlnitRsp 
©inline SpSat_blscaie.e SpSatUpdateRsp 
©inline SpSat_blscale.e SpSatON 
©inline Prescan.e PScore 

45 

/• -PSDINIT- •/ 
int 
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psdinitQ 
{ 

int simod_acqs; 
im i; 

5 

f* Initialize everything to a known state •/ 

setrfconfig((shon)5); /* only activate rhol V 
setssitime(time_ssi / GRAD_UPDATE_TIME); 
rspqueueinit(queue_size); /* was 200 •/ 

10 scopeon(&seqcore); /• Activate scope for core •/ 

syncon(&seqcore); /* Activate sync for core */ 
syncoff(&seqpass); /* Deactivate sync during pass •/ 
seqCount = 0; . /* Set SPGR sequence counter V 
settriggerarray((short)slquanth rsptrigger); 

15 setrotatearray((shon)slquantI, rsprot); 

chopamp = ia_rfl; 
setiamp(chopamp. &rfL 0); 
for (i = 0; i < opetl: i— r) { 

setrfltrsfscanslot. &echolfil); 

20 } 

/* sat / chemsat */ 

SpSatInitRsp(vrgsat. sp_satindex. sp_satcardJoc. 0); 
SpSat_set_satl_matrix(rsprot_orig, rsprot. opslquant • opphases, 
25 sat_rot_matrices. sat_rot_ex_num. sat_rot_df_num. 

sp_satcard_Ioc. 0); 
setrotatearray((SHORT)(opsIquant * opphases). rsprot); 
CsSatMod(cs_satindex); 

30 /* setup pass arrays */ 

slmod_acqs = opslquant % acqs: 
for (i = 0: i < acqs: i-r-) { 

sic _in_acq[i] = opslquant / acqs; 
if (i < s!mod_acqs) 
35 slcjn_acq[i] += I; 

if 0 — 6) 

ac(Lptr[i] = 0; 
} else { 

acqj)tr[i] = acqjptr[i - 1] 
+ slc_in_acq[i - 1]; 

} 

} 



40 



/• asym FOV - overrides previous code •/ 
45 /* set up phase offset arrays •/ 

for (i = 0; i < opslquant * opfphases; H-+) { 
if (rsp_shift[i].rspphasoff >= 0) { 
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phase_off[i].ysign = -1; 
} else { 

phase_off[i].ysign = 1; 

} 

5 /* phase offset increment •/ 

yoffsl = 0.5 + fabs(FS_2PI • rsp_shift[i].rspphasoff / 
(opfov a opphasefov*nop)); 
/• offset in ranee */ 

phase_off[i].yoffs = (yoffsl + FS_2PI + FS PI) % FS 2PI 
10 - FS PI; 

} 

sp_sat index = 0; 
if (opsat = PSDON) { 
15 if (ssivector(ssisat) = FAILURE) return FAILURE; 

} else { 

if (ssivector(dummyssi) = FAILURE) return FAILURE; 

} 

' 20 return SUCCESS: 

) 

int 

mps2() 
25 { 

rspent = L_MPS2; 

rspdda = 2; 

rspbas = 0; 

rspvus = 30000; 
30 rspgyl = 0: 

rspnex = 2; 

rspesl = pre_slice: 

rspasl = pre_slice; 

rspslq = 1; 
35 rspsct = 0; 

cstun = 1 ; /• turn chemsat on */ 

psdinit(); 
SpSat_Saton(0); 
40 setperiod(optr - seqtime, &seqcore. 0); 

boffset(off_seqcore); 
pscore(); 

return rspexit(); 
45 } /• mps2() •/ 
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10 



15 



20 



25 



30 



35 



40 



aps2() 
i 

rspent = L_APS2; 
rspdda = 2; 
rspbas = 0; 
rspvus = 1024; 
rspgyl = 0; 
rspnex = 2; 
rspesl = pre_slice; 
rspasl = pre_slice; 
rspslq = I; 
rspsct = 0; 

cstun =1; /• turn chemsat on •/ 

psdiniU); 
SpSat_Saton(0); 

sctperiod(optr - seqtime, &seqcore. 0); 

boffset(off_seqcorc): 

pscore(); 

return rspexit(); 
} /• aps2() V 

/• -SCAN- •/ 
int 

scan() 
{ 

short pause; 
int pass; 

rspent - L_SCAN; 
rspsct = 0; 
rspslq = slquantl; 

cstun =1; /• turn chemsat on */ 

psdinit(); 
SpSat_Saton(0); 

setperiod(deadtime. &seqcore, 0); 

for (pass = 0; pass < acqs; pass++) { 
/• scan •/ 

boffset(off_seqcore); 

passjp(pass); 

/• pass packet •/ 

boffsetfoff_seqpass); 

if (pass = (acqs - 1)) { /• last pass, eoscan •/ 
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45 



setwamp(SSPD + DABPASS + DABSCAN. &pass_pulse. 2); 
} else { /* eopass •/ 

setwamp(SSPD + DABPASS, &pass_pulse. 2); 

} 

/* pause •/ 

if (pass = (acqs - 1)) { 

pause = MAYJ>AUSE; 
} else { 

if (((pass + 1) % slicecnt) = 0) { 

pause = MUST_PAUSE; /• pause if desired •/ 
} else { 

pause = AUTO J>AUSE: f or if required •/ 

} " 

} 

$p_sat_index = 0; 
stanseq(0, pause); 



} 



} 

return rspexiif); 



int 

rot_corr(sl) /• rotate •/ 
int si; 



{ 



double th. cs. sn; 
float view_rot[3][3]; 
shon net_rot[9]; 
int i. j, k. ix: 

view_rot[0][0] = rll: view_rot[0][l] = rl2; viewrot[0][2] = r!3; 

view_rot[I][0] = r21: view_roi[l][l] = r22; view_rot[l][2] = r23; 

view_rot[2][0] » r31; view_rot[2][l] = r32; view_rot[2][2] = r33; 

for (ix = 0; ix < 9; ix-n-) net_rot[ix] = 0; 
for (i * 0; i < 3; i-r+) { 

for(j = 0:j<3:j+-) { 
ix = j + 3 • i; 
for (k = 0; k < 3; k++) { 

net_rot[ix] += (shon) 
/• (orig_rot[i][k] • view_rot[k][j]); */ 

(view_rot[i][k] • orig_rot[k][j]); 

} 

} 

} 

setrotate(net_rot, si); 



return 0; 
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/• -CORE- •/ 
5 int 

pass_lp(pass) /* spin echo core ... multi-pass */ 

int pass; 

{ 

int view; /* view index */ 

10 int shot; /• excitation index •/ 

int excitation; /* NEX index •/ 

int i; 

for (shot = -dda: shot < nshot: shot++) { 
15 for (i - 0; i < opetl; i^) { 

if (dbg = 0) { 

view - shot; 
} else { 

view = i V nshot - shot: 

20 } 

/* phase encoding view */ 
if (shot < 0) view = 0; /* if DDA •/ 
setiamp(viewtable[view + 1], &gyl, i); 
setiamp(viewtable[view + 1]. &evla. i); 

25 } 

/* excitation loop */ 

for (excitation = 0; excitation < opnex: excitation^) { 
/* DABOP and chop V 
30 if (excitation = 0) { 

dabop = 0; 
} else { 

dabop = 1 + 2 * (excitation % 2); 

} 

35 getiampO&chopamp, &rfl, 0); 

setiamp(-chopamp. &rfl, 0); 



40 



/* scan •/ 
slicejp(pass, shot); 



} 



slicejp(pass. shot) 
45 int pass, shot; 
{ 

int view; /* view index V 
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float phs; 

int yres_phase; 

int echo, \J* echo index */ 

int sK slice: /• slice index ■/ 

int dummy_flag. dda_flag; 

if (shot >= 0) { 

dda_flae « PSDJ3N; 
}e!se { 

ddajlae = PSD_OFF; 

} 

dummy_flag = PSDJDN; 

for (si = 0; si < rspslq; sl++) { 

/* rotate •/ 

if (rot_fix) rot_corr(sl); 

/* scope trigger V 

if (rspsct = -1 || si = rspsci) { 

scopeonf &seqcore): 
} else { 

scopeofR&seqcore): 

} 

/* dummy slice for multi-pass */ 
if (si >~ slc_in_acq[pass]) { 

dummy_flag = PSDJ3FF; 

seueos((short)EOS_DEAD. &rho tdO. 0); 
} else { 

dummy Jlag - PSDJDN; 

setieos((shon)EOS_PLAY, &rho_tdO, 0); 

} 

/• set freq and go •/ 

slice = acq_ptr[pass] + si; 

if (slice >= opslquant) slice ~ 0; 

/• RF freq and phase •/ 

setfrequency(rfl_freq[slice], &rfl, 0); 
setfrequency(r£2_freq[slice], &rfZ 0); 
setfrequency(rf2_freq[siice], &r£22, -1); 
setphase(rfl_phase[slice], &rfl, 0); 
setphase(rf2_phase[slice], <&rf2, 0); 
setphase(rf2_phase[slice], &rf22, -1); 

/• receiver freq and phase. DAB packet (DAB view etc) */ 
if (dda_flag && dummy_flag) { 
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acq_echol = (int)DABON; 
} else { 

acq_echol = (int)DABOFF; 

} 

for (i = 0; i < opetl: i^) { 
if (dbg = 0) { 

echo = i; /• CPMG */ 

view = shot; 
} else { 

echo = 0; /• RARE •/ 

view = i • nshot + shot; 

} 

loaddab(&echol[i], si, echo, dabop. 

view + 1, acq_echol); 
yres_phase = phasejjfljslicej.ysign 

♦ ((view * phase offjslicej.yoffs 
+ 3L * FS_PI) 
% FS_2PI 
- FS_PI); 
setiphase(yres_phase, &echol[i], i); 
setfrequency(receivejreql[slice]. &echol[i], i); 



/• stanseq() •/ 

sp_sat_index = si; 

stanseq(sK (short)MAY_PAUSE); 

syncoff(&seqcore); 

J 

} 

/• -PSCORE- •/ 
int 

pscore() 
{ 

int view, i; 

setfrequency(rfl_freq[rspesl], Arfl. 0); 
setfrequency(rG_freq[rspesl], &rf2, 0); 
setfrequencyCreceive^freqllrspasl], &echoi[0], 0); 

dabop = 0; 

for (i = 0; i < opetl; i-w-) { 

loaddab(&echoI[i], 0, 0, dabop, 0, DABOFF); 

} 

setiamp(0, &gyl, 0); 

for (view = 0; view < rspdda: view-M-) { 
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sp_sat_index = 0; 
startseq(0, (shon)MAY_PAUSE); 
gctiamp(&chopamp. &rfl, 0); 
setiamp(-chopamp. &rfh 0); 

} 

for (view = 0; view < rspvus; view++) { 

Ioaddab(&echol[0], 0, 0, dabop. view + 1, DABON); 

sp_sat_index = 0; 

stanseq(0, (shon)MAY_PAUSE); 

getiamp(&chopamp, &rfl % 0); 

setiamp(-chopamp, &rfl, 0); 

} 

} /* pscore() •/ 



dsl/206215.2 
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35 

WHAT IS CLAIMED IS: 

1. A method of obtaining imaging scans of a subject 
taken over different time periods comprising: 

(a) taking at a first time period, using an imager, a 
5 baseline imaging scan that is relatable to a personal 

coordinate system (PCS) and a first machine coordinate system 
(MCS) , which PCS can be regenerated from the subject in a 
substantially reproducible manner ; 

(b) taking at a different time period, using a second 
10 imager, at least one follow-up imaging scan, which scan is 

relatable to the PCS and a second MCS; 

(c) manipulating the second MCS, such that the 
relationship between the second MCS and the PCS is 
substantially the same as the relationship between the first 

15 MCS and the PCS, to permit the taking of at least one 
additional follow-up imaging scan that can be superimposed on 
the baseline imaging scan. 



2. The method of claim 1 which further comprises taking 
20 one or more additional follow-up imaging scans that can be 

superimposed on the baseline imaging scan. 

3. The method of claim 1 in which the second imager of 
step (b) is the same as the imager of step (a) . 

25 

4. The method of claim 1 in which the second imager of 
step (b) is different from the imager of step (a) . 

5 . The method of claim 1 in which the MSC of the imager 
30 of step (a) is fixed. 



6. A method of obtaining imaging scans of a subject 
comprising : 

(a) providing a subject with a non-invasive stereotactic 
35 device that is positioned reproducibly on a subject and which 
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establishes a personal coordinate system (PCS) associated with 
the subject which is independent of a machine coordinate system 
(MSC) associated with an imager; 

(b) taking, using an imager having an MCS, an imaging 
5 scan of the subject including the stereotactic device to 

establish the PCS of the subject; 

(c) manipulating the MCS of the imager to bring the MCS 
in substantial alignment with the PCS of the subject ; 

(d) taking one or more additional imaging scans of the 
subject with the MCS of the imager substantially aligned with 
the PCS of the subject, to obtain a first set of imaging scans. 

7. The. method of claim 1 which further comprises 
repeating steps (a) to (d) at a second time period, using a 
second imager, to obtain a second set of imaging scans. 

8. The method of claim 7 which further comprises 
comparing at least one imaging scan of the first set with at 
least one imaging scan of the second set. 

9. The method of claim 8 which further comprises noting 
any previously undetected anatomical feature of the subject. 

10. The method of claim 8 which further comprises 
observing any changes in any previously detected anatomical 
feature of the subject. 

11. The method of claim 9 in which said anatomical 
feature comprises a tumor. 

12. The method of claim 10 in which said anatomical 
feature comprises a tumor. 

13 . The method of claim 7 in which the imagers are the 

same . 
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14 . The method of claim 7 in which the imagers are 
different . 

5 15. The method of claim 6 in which the imager is a 

magnetic resonance imager. 

16. The method of claim 6 in which the stereotactic 
device is positioned reproducibly on the patient's head. 

10 

17. The method of claim 7 in which the second time period 
represents an elapsed time from the taking of the first set of 
imaging scans to the taking of the second set of imaging scans 
of about one day to about one week, about one week to about one 

15 month, about one month to about six months and about six months 
to about one year. 

18. The method of claim 7 in which the second time period 
represents an elapsed time from the taking of the first set of 

20 imaging scans to the taking of the second set of imaging scans 
of about one year to about five years. 

19. A method of obtaining imaging scans of a subject 
comprising : 

25 (a) providing a subject with a non-invasive stereotactic 

devi ce that is positioned reproducibly on a subject and which 
establishes a personal coordinate system (PCS) associated with 
the subject which is independent of a machine coordinate system 
(MSG) associated with an imager; 

30 (b) taking, using a first imager having a first MCS, at 

least one imaging scan of the subject including the 
stereotactic device to establish the PCS of the subject and to 
relate the PCS of the subject to the first MCS of the first 
imager; 
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(c) taking, using a second imager having a second MCS, 
at least one imaging scan of the subject including the 
stereotactic device to reestablish the PCS of the subject and 
to relate the PCS of the subject to the second MCS of the 

5 second imager; 

(d) manipulating the second MCS , such that the PCS is 
related to the second MCS in substantially the same way as the 
PCS is related to the first MCS. 

10 20. The method of claim 19 which further comprises taking 

one or more additional imaging scans of the subject with the 
second imager. 

21. The method of claim 19 which includes substantially 
15 aligning the second MCS with the first MCS. 

22. The method of claim 19 in which the imagers are 
different / 

20 23. The method of claim 22 in which the first imager is 

a CT scanner. 

24. The method of claim 22 in which the second imager is 
an MR imager. 

25 

25. The method of claim 22 which further comprises 
forming a composite image including information from at least 
one imaging scan taken using the first imager and information 
from at least one imaging scan taken using the second imager. 

30 

26. A stereotactic device for use with an imager 
comprising a frame equipped with localizing means and affixing 
means , 

the localizing means comprising one or more localizing 
35 arrays that provide one or more imager detectable signals, 
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the affixing means comprising non- invasive fittings for 
placement about the periphery of the subject and which permit 
the reproducible positioning of the frame on the subject, 

from which signals is derived a personal coordinate system 
5 that serves as a reference coordinate system for imaging scans 
taken of a subject on which the frame is positioned. 

27. The device of claim 26 which facilitates the 
stereotactic examination of the head of the subject. 

10 

28. The device of claim 26 in which the localizing array 
comprises one or more reference elements. 

29. The device of claim 28 in which the reference 
15 elements are radio-opaque. 

30. The device of claim 29 in which the reference 
elements are detectable by an MR imager. 

20 31. The device of claim 26 in which the localizing means 

comprises four or more localizing arrays each comprising a 
reference element , 

each reference element defining a point in space, three 
of which points define a unique plane and the fourth point 

25 lying outside the unique plane,. 

32. The device of claim 26 in which the localizing means 
comprises three localizing arrays each comprising a pair of 
reference elements . 

30 

33. The device of claim 32 in which the paired reference 
elements are arranged in a spaced-apart W X" configuration. 

34. The device of claim 33 in which the paired reference 
3 5 elements are orthogonal to one another. 
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35. A stereotactic device for use with an imager 
comprising a frame equipped with localizing means and affixing 
means , 

5 the localizing means comprising three or more localizing 

arrays each equipped with two or more reference elements that 
together provide six or more imager detectable signals, 

the affixing means comprising three or more non- invasive 
fittings for placement on the subject's nasion and in or about 
10 the subject's ears and which permit the reproducible 
positioning of the frame on the subject, 

from which signals is derived a personal coordinate system 
that serves as a reference coordinate system for imaging scans 
taken of a subject on which the frame is positioned. 

15 

36. The device of claim 35 in which the reference 
elements comprise elongate components. 

37. The device of claim 36 in which two elongate 
20 components are arranged in an M X" configuration and in which 

the two elongate components are related to one another by a 
predetermined angle. 

38. The device of claim 37 in which the predetermined 
25 angle is about 90E. 

39. The device of claim 35 in which a localizing array 
is disposed proximate to an ear fitting. 

30 40. The device of claim 35 in which the reference 

elements comprise radiographically opaque or semi -opaque 
material . 

41. The device of claim 40 in which the material 
35 comprises steel or doped water. 
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42. The device of claim 35 which further comprises a 
securing means that facilitates the securing of the frame to 
the subject. 

5 

43. The device of claim 42 in which the securing means 
comprises an inelastic or elastic component and facilitates the 
securing of the frame to the head of the subject. 

10 44. The' device of claim 26 in which the one or more 

localizing arrays provide one or more signals that are 
detectable by MR, CT f PET, SPECT, EEG, or MEG . 

45. The device of claim 26 in which the imager detectable 
15 signals are provided by a digitizer used in conjunction with 

the one or more localizing arrays. 

46. A method of spatially aligning at least two 
radiographic scan images of the head of a subject taken by a 

20 radiographic scanning device comprising: 

(a) radiographically scanning the head of a subject 
wearing the stereotactic device of claim 1 to provide a first 
radiographic scan image containing a plurality of first 
reference points defining a first personal coordinate system; 

25 (b) relating said first personal coordinate system with 

a personal coordinate system obtained from a prior radiographic 
scan of said subject reproducibly wearing said stereotactic 
device; 

(c) adjusting said radiographic scanning device to align 
30 said first and prior personal coordinate systems; 

(d) obtaining one or more additional radiographic scans 
of said subject with the adjusted radiographic scanning device. 
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